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ABSTRACT 


The  object  of  the  thesis  is  development  of  a  method  for 
describing  the  behavior  of  a  ship's  material  condition  over 
time  so  that  the  way  in  which  various  maintenance  strategies 
affect  this  material  condition  may  be  examined.   The  thesis 
demonstrates  the  manner  in  which  these  strategies  may  then  be 
evaluated  based  on  total  expected  cost  of  a  fleet  of  ships. 

Ship  failure  and  deterioration  behavior  is  modeled  as  a 
transient  Markov  process,  with  a  set  of  increasingly  complex 
models.   The  most  comprehensive  of  these  allows  such  behavior 
to  depend  on  the  ship's  present  condition,  its  age,  and  the 
time  since  it  was  last  overhauled.   In  the  models  a  ship  is 
treated  as  a  single  component,  rather  than  an  aggregation  of 
several  smaller  components,  and  it  is  assumed  that  the  ship's 
material  condition  can  be  categorized  in  a  discrete  manner. 

Use  of  the  models  is  demonstrated  by  the  use  of  two 
example  cases,  one  of  which  is  based  on  actual  maintenance 
data  for  Navy  ships . 
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I .   INTRODUCTION 

A.   OBJECTIVES: 

It  is  the  intention  of  this  thesis  to  develop  an  approach 
for  examining  the  behavior  of  a  ship's  material  condition  over 
time  and  the  manner  in  which  that  behavior  is  affected  by  var- 
ious alternative  maintenance  policies.   We  will  do  this  under 
the  premise  that  a  trade-off  exists  between  the  investment 
cost  of  a  fleet  of  ships  and  maintenance  expenditures  for  the 
fleet  under  various  maintenance  policies,  given  an  operational 
requirement  of  the  form  that  some  number  of  ships  x  must  be 
available  for  operations  with  probability  p.   For  a  specified 
x  and  p,  our  goal  will  be  to  determine  that  combination  of 
ship  acquisition  and  overhaul  policies  with  the  lowest  expec- 
ted cost,  which  obeys  this  constraint. 


B .   BACKGROUND : 

In  recent  years  nearly  20  percent,  or  approximately  a 
billion  dollars  a  year,  of  the  Navy's  entire  operation  and 
maintenance  budget  has  been  spent  on  the  maintenance  of  the 
Navy's  ships.1   This  does  not  include  the  cost  of  naval  man- 
power, the  cost  of  support  ships  and  advance  support  bases, 
nor  the  capital  cost  of  naval  shipyards,  all  of  which  contri- 
bute significantly  to  the  real  cost  of  ship  maintenance.   Nor 
does  it  include  any  additional  investment  in  ships  that  may 
result  from  choosing  one  maintenance  policy  over  another. 

A  significant  part  of  the  cost  of  ship  maintenance  is 
incurred  during  the  ship's  regular  overhaul.   A  regular  over- 
haul is  a  scheduled  period  (under  present  Navy  policy) ,  occur- 
ring at  specified  intervals,  during  which  a  ship  is  assigned 
to  a  shipyard  for  work  of  a  relatively  major  nature.   It 
includes  maintenance,  both  corrective  and  preventive,  as  well 
as  (usually)  some  items  of  modernization. 

Because  a  significant  portion  of  the  maintenance  dollar 
is  spent  on  regularly  scheduled  ship  overhauls,  a  great  deal 
of  attention  has  been  focused  on  improving  the  effectiveness 
of  overhauls,  i.e.,  attempting  to  get  maximum  benefit  from  the 


*See  U.S.  Office  of  Management  and  Budget,  Budget  of  the 
United  States  Government  and  appendices  thereto,  1969-1972, 
Government  Printing  Office,  Washington,  DC;  and  Department  of 
the  Navy,  Naval  Ship  Systems  Command,  The  Cost  Estimate  Fact 
Book,  Volume  I:      Systems  and  Methods,  Washington,  DC,  1969. 


10 
work  accomplished  relative  to  the  amount  of  time  and  money 
invested.   This  effort  has  generally  proceeded  along  two 
lines:   (1)   By  attempting  to  shorten  overhauls  while  main- 
taining the  same  amount  of  work  accomplished  by  improved 
planning,  skill  utilization,  material  procurement,  etc.,  by 
such  efforts  as  the  Compressed  Regular  Overhaul  (CROH)  Pro- 
gram;  (2)   By  examining  the  appropriateness  of  the  length 
between  regular  overhauls  and  the  manner  in  which  a  ship's 
material  condition  deteriorates  between  regular  overhauls. 
It  is  in  this  latter  area,  perhaps,  that  the  most  effort  has 
been  exerted. 

In  this  area,  again,  two  basic  approaches  seem  to  have 
been  taken:   The  first  of  these  is  to  examine  separately  the 
failure  behavior  of  the  various  components  which  make  up  a 
ship  and  then  attempt  to  synthesize  this  into  a  description  of 
the  behavior  of  the  material  condition  of  the  overall  ship. 
To  conduct  such  an  analysis  for  just  one  component,  assuming 
the  data  to  support  the  analysis  is  available,  is  a  rather 
extensive  task.   Frequently,  sufficiently  complete  or  accurate 
data  to  adequately  support  such  an  analysis  is  not  available 
for  the  component.   In  considering  an  entire  ship  made  up  of 
several  hundred  subsystems  the  magnitude  of  the  job  becomes 
quite  formidable.   In  the  Navy,  this  is  compounded  by  the  lack 
of  commonality  between  components  of  ships  of  the  same  type, 
or  even  of  the  same  class  of  ship  within  a  type.   This  has 
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resulted  in  part  from  the  Navy's  policy  (in  the  past,  at 
least)  of  letting  contracts  for  ships  of  a  given  class  to  sev- 
eral different  building  yards  without  requiring  commonality  of 
equipment. 

A  more  promising  path  to  take  (and  the  second  basic 
approach  to  the  problem),  at  least  in  the  short  run,  would 
seem  to  be  to  treat  the  whole  ship  as  a  single  component. 
Indeed,  using  this  approach,  Farrar  and  Apple2  were  reason- 
ably successful  in  developing  useful  formulas  for  repair 
requirements  as  well  as  information  on  the  difference  between 
scheduled  and  unscheduled  repairs.   This  was  done  in  spite  of 
much  of  their  data  being  incomplete  and  somewhat  inaccurate. 
Their  work  tends  to  support  the  basic  feasibility  of  this 
second  approach,  in  that  a  great  deal  of  useful,  yet  relative- 
ly inexpensive  information  can  be  derived  by  considering  the 
ship  as  a  single  unit. 


2See  Donald  E.  Farrar  and  Robert  E.  Apple,  "Some  Factors 
That  Affect  the  Overhaul  Cost  of  Ships:   An  Exercise  in  Sta- 
tistical Cost  Analysis,"  Naval  Research  Logistics  Quarterly, 
vol.  10,  no.  4,  December  1963,  pp.  335-368,  and  "Economic  Con- 
siderations in  Establishing  an  Overhaul  Cycle  for  Ships:   An 
Empirical  Analysis,"  Institute  of  Naval  Studies  Research  Con- 
tribution No.  7,  April  1964,  Naval  Engineer' s  Journal,  vol. 
77,  no.  1,  February  1965,  pp.  69-78. 
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C.   THE  APPROACH: 

Here  we  pursue  the  "whole  ship"  approach,  and  develop  an 
increasingly  general  set  of  stochastic  models  describing  the 
behavior  of  a  ship's  overall  material  condition  as  the  ship 
ages  and  as  alternate  overhaul  strategies  are  employed.   These 
models  will  be  designed  to  be  intuitively  appealing  to,  and 
in  fact  to  make  use  of,  the  informed  judgment  of  Navy  decision 
makers.   A  model  which  a  decision  maker  does  not  understand, 
or  which  he  finds  at  odds  with  his  intuition,  is  unlikely  to 
see  much  application. 

Once  a  model  has  been  developed  through  which  a  reason- 
able approximation  of  the  behavior  of  a  ship's  material  con- 
dition over  its  life  can  be  made,  a  number  of  different 
investment  and  maintenance  policies  will  be  studied  for  their 
effect  on  the  total  cost  of  meeting  a  specified  operational 
requirement.   In  particular,  once  the  behavior  of  a  single 
ship  has  been  modeled,  information  can  then  be  determined  for 
an  entire  fleet  of  such  ships.   For  example,  the  expected  num- 
ber of  ships  in  overhaul  at  a  given  time  for  a  given  fleet 
size  or,  given  the  requirement  for  maintaining  a  specified 
number  of  operational  ships  with  a  specified  probability,  the 
required  fleet  size  can  be  determined.   We  will  show  this  in 
the  thesis. 

In  this  thesis  we  will  consider  cost  as  belonging  to  two 
categories:   overhaul  cost,  and  the  ship  acquisition  and  life 
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cycle  costs  not  including  overhaul  cost  (reduced  to  present 
value) .   The  latter  will  henceforth  be  referred  to  as  "ship 
cost".   In  order  to  be  valid  for  decision  making,  this  must 
include  the  present  value  of  the  ship's  entire  life  cycle 
costs  less  overhaul  costs,  including  cost  of  crew,  etc.,  as 
well  as  the  actual  acquisition  cost.   If  we  know  or  can  reas- 
onably estimate  ship  cost  as  well  as  overhaul  cost,  and  if  we 
can  determine  the  expected  number  of  overhauls  the  ship  will 
experience  under  a  given  policy,  then  we  can  determine  total 
expected  lifetime  cost. 

The  models  developed  in  the  thesis  will  include  the  capa- 
bility to  investigate  the  effect  on  total  fleet  cost  of  var- 
ious overhaul  policies  and  the  required  fleet  size  (number  of 
ships)  resulting  from  these  policies. 

We  will  model  ship  failure  and  deterioration  behavior  as 
a  transient  Markov  process.   Beginning  with  an  extremely  sim- 
ple process  we  will  elaborate  it  until  it  has  the  ability  to 
reflect  observed  behavior  in  a  reasonably  complete  manner. 
Once  this  has  been  done,  data  extracted  from  actual  ship  main- 
tenance records  will  be  introduced  and  use  of  the  most  com- 
prehensive model  as  an  aid  to  decision  making  will  be  demon- 
strated. 

This  approach  to  modeling  a  ship's  change  in  material 
condition  over  time  represents,  it  is  felt,  a  somewhat  unique 
conceptualization  of  the  problem.   Although  Farrar  and  Apple 
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treated  the  ship  as  a  single  component,  their  approach  depen- 
ded primarily  on  the  use  of  regression  analysis. 

The  model  should  be  applicable  to  both  naval  and  commer- 
cial ships,  at  least  in  basic  concept.   However,  in  the 
thesis,  where  specifics  are  required  in  either  developing  or 
demonstrating  the  model,  they  will  be  derived  from  and  geared 
to  destroyer-type  ships  of  the  U.S.  Navy. 
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D.   SUMMARY: 

The  remaining  chapters  of  the  thesis  are  summarized  as 
follows:   In  Chapter  II  we  first  establish  the  underlying 
assumptions  for  the  model.   We  then  develop  a  simple  model  in 
which  two  categories  of  ship  material  condition  are  assumed 
(consider  these  as  "failed"  and  "not  failed")  and  in  which 
the  probability  of  failure  is  dependent  solely  on  the  ship's 
age.   Chapter  II  concludes  with  a  hypothetical  example  demon- 
strating use  of  the  model  as  developed  at  that  point. 

In  Chapter  III  we  introduce  into  the  model  the  assumption 
that  not  only  will  the  probability  of  failure  change  with 
ship  age,  but  with  the  length  of  time  since  the  ship  was  last 
overhauled  as  well.   We  then  expand  our  binary  failure  assump- 
tion (two  categories  of  material  condition)  to  one  which 
allows  for  four  categories  of  ship  material  condition,  in 
which  the  worst  of  these,  condition  4,  corresponds  to  failure 
in  the  binary  model.   We  conclude  Chapter  III  with  an  example, 
using  our  four  condition  model,  based  on  data  extracted  from 
the  maintenance  records  of  the  Commander,  Cruiser-Destroyer 
Force,  Atlantic  (COMCRUDESLANT) . 

Chapters  II  and  III  of  necessity  deal  to  a  certain  extent 
with  the  mechanics  of  computation  for  the  model.   It  is  not 
until  we  reach  Chapter  IV  that  we  are  actually  able  to  apply 
the  model  to  a  comparison  of  alternative  overhaul  policies, 
utilizing  the  example  of  Chapter  III  as  a  basis.   In  Chapter 
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IV  we  examine  the  effect  of  overhaul  policy  on  the  number  of 
ships  necessary  to  meet  a  given  operational  requirement  with 
a  specified  probability.   Chapter  IV  compares  overhaul  poli- 
cies of  several  specified  intervals  with  the  policy  of  over- 
hauling a  ship  only  upon  failure,  i.e.  a  "demand"  dependent 
policy.   Chapter  IV  further  considers  the  effect  of  varying 
the  specified  operational  requirement  and  the  effect  of  im- 
puting some  penalty  cost  to  an  unscheduled  overhaul,  i.e., 
one  which  occurs  on  demand  due  to  a  ship  failure. 

In  Chapter  V,  after  discussing  some  possible  directions 
for  further  work  in  this  area,  we  arrive  at  the  following 
conclusions  resulting  from  the  thesis: 

(1)  That  the  behavior  of  a  ship's  material  condition 
over  time  can  be  successfully  modeled  by  a  transient  Markov 
process . 

(2)  That  such  a  model  can  be  of  use  to  Navy  decision 
makers  in  comparing  alternative  overhaul  policies. 

(3)  That  actual  implementation  of  the  model  will  require 
initially  a  major  data  collection  and  reduction  effort. 
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II.   DEVELOPMENT  OF  THE  BASIC  MODEL 

A.   ASSUMPTIONS: 

The  basic  assumptions  underlying  the  models  used  in  this 
report  are:   (1)   that  the  change  in  a  ship's  material  condi- 
tion over  time  may  be  described  by  a  Markov  process,  (2)   that 
the  ship  can  be  classified  into  one  of  two  or  more  discrete 
categories  of  ship  material  condition,  and  (3)   that  some 
reasonable  approximation  of  the  probabilities  of  remaining  in 
the  same  or  changing  to  a  different  material  condition  cate- 
gory as  the  ship  ages  can  be  determined. 

A  first  order  Markov  process  is  one  in  which  only  the 
last  state  occupied  by  the  process  is  relevant  in  determining 
its  future  behavior.   Under  the  Markovian  assumption  the  prob- 
ability of  making  a  transition  to  each  state  of  the  process 
depends  only  on  the  state  presently  occupied. 3 

As  noted  earlier,  a  ship,  especially  a  warship,  is  suf- 
ficiently complex  that  its  overall  material  condition  is  com- 
prised of  the  net  effect  of  the  material  condition  of  a  large 
number  of  component  systems  and  subsystems.   Further,  it  can 
be  significantly  affected  by  a  number  of  diverse  factors,  such 
as  the  intensity  of  operations,  the  funds  available  for  re- 
pairs and  spare  parts,  the  qualification  and  motivation  of  the 


'See  Ronald  A.  Howard,  Dynamic  Probabilistic  Systems,  Vol- 
ume I_:   Markov  Models,  John  Wiley  and  Sons,  New  York,  1971, p.3. 
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operating  crew,  etc.   Consequently/  a  means  of  precisely 
determining  a  ship's  overall  material  condition  has  thus  far 
been  somewhat  elusive.   Efforts  to  devise  some  sort  of 
material  condition  index  have  been  initiated  and  are  contin- 
uing; however,  no  such  usable  index  exists  at  this  time. 

It  seems  rather  that  the  best  means  to  try  to  determine 
a  ship's  overall  material  condition  is  to  ask  the  operator/ 
maintainer,  in  the  present  case  the  type  commander.   Such  an 
approach  is  necessarily  subjective  and  is  more  likely  to  be 
qualitative  than  quantitative  in  nature.   The  type  commander, 
supported  by  his  staff,  is  nonetheless  likely  to  have  a  good 
"feel"  for  the  condition  of  his  ships  since  he  is  in  frequent 
contact  with  them  and  has  the  benefit  of  experience  with  var- 
ious ships.   Based  on  this,  the  model  assumes  that  a  type  com- 
mander, utilizing  the  collective  knowledge  and  experience  of 
his  operating  and  maintenance  staffs,  can  (and  in  practice 
does,  although  perhaps  not  formally)  categorize  his  ships  into 
two,  three  or  four  different  levels  of  material  condition. 
Some  artificiality  exists  in  that  the  model  assumes  discrete- 
ness where  in  fact  a  continuum  exists.   Nonetheless,  exper- 
ience has  indicated  that  a  type  commander  can  fairly  well  tell 
which  of  his  ships  he  considers  "good",  which  are  "bad",  and 
which  he  would  group  in  the  middle  as  not  belonging  to  either 
extreme  category. 

Initial  estimates  of  the  probability  of  transition  from 
one  of  these  categories  to  another,  as  the  ship  ages,  must  be 
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made.   Obviously  the  best  means  of  doing  this  is  to  observe 
the  behavior  of  a  large  number  of  similar  ships  over  a  period 
of  time.   Unfortunately,  in  the  Navy  to  date  there  has  been 
little  or  no  formal  categorization  of  total  ship  material  con- 
dition as  described  above,  so  no  past  records  are  available 
from  which  estimators  of  transition  probabilities  may  be 
determined  directly,  although  a  similar  system  is  in  effect  in 
the  Navy  to  denote  the  seriousness  to  the  ship's  operating 
capability  of  material  casualties  to  individual  components 
(part  of  the  Casualty  Reporting,  or  CASREPT,  system) .   Conse- 
quently, the  best  approach  would  seem  to  be  Bayesian  in 
nature,  where  prior  estimates  of  the  behavior  of  a  ship's 
material  condition  are  established  based  on  the  best  data  and 
informed  opinion  available,  and  then  are  updated  as  better 
information  becomes  available.   A  discussion  of  how  such  prior 
estimates  might  be  made  will  appear  later. 


20 

B.   THE  TWO  CONDITION,  AGE  DEPENDENT  MODEL: 

Given  the  foregoing  assumptions,  we  will  begin  our  ana- 
lysis by  considering  a  relatively  simple  case.   In  this  sim- 
plest case,  two  possible  conditions  are  assumed:   these  can 
be  considered  "good"  and  "bad",  "operating"  and  "failed",  etc. 
They  will  be  designated  as  "1"  and  "2",  respectively. 

Intuitively,  it  seems  reasonable  to  expect  a  higher  prob- 
ability of  failure  for  an  old  ship  than  for  a  new  one.   While 
this  may  or  may  not  be  true,  it  would  be  desirable  if  the 
model  were  able  to  accomodate  such  a  possibility.   In  order 
adequately  to  describe  a  state,  therefore,  both  the  ship's 
condition  and  age  must  be  specified.   Let  the  state  indices  be 
age  t  and  condition  k;  the  state  then  is  described  by  the 
ordered  pair: 

(k,t)  =  ship  is  t  units  old  and  in  condition  k, 
keK,  teH, 
where  K  =  {1,2}  for  a  two  condition  model, 

and  H  =  {1,2, ...,T},  T  being  the  life  (i.e.,  maximum 
age  the  ship  will  reach)  of  the  ship  in 
terms  of  the  basic  time  unit. 
The  assumption  is  made  that  the  ship  will  occupy  the 
state  for  all  of  a  given  time  increment.   This  should  be  reas- 
onably valid  if  the  time  increment  is  relatively  short.   In 
fact,  it  could  correspond  to  a  periodic  appraisal  of  the  con- 
dition of  a  fleet  of  ships,  so  that  transition  would  take 
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place  only  when  the  ships'  condition  was  reappraised.   In  this 
simple  model,  it  is  assumed  that  the  decision  maker,  upon 
finding  a  ship  in  condition  2,  would  immediately  take  correc- 
tive action  consisting  of  an  "overhaul".   In  other  words,  a 
"demand"  (for  overhaul)  condition  would  be  said  to  exist  any 
time  the  ship  reached  condition  2.   The  objective  of  over- 
hauling a  ship  presumably  is  to  improve  its  material  condi- 
tion.  It  will  therefore  be  assumed  that  the  overhaul  will 
with  certainty  return  the  ship  to  condition  1.  Her  AtW 

Since  an  overhaul,  under  the  latest  assumption,  in  effect 
represents  a  transition  from  condition  2  to  condition  1,  it  is 
convenient  to  take  as  the  basic  time  unit  a  period  equal  to 
the  approximate  length  of  overhaul  of  a  ship,  or  a  period  of 
which  length  of  overhaul  is  a  multiple.   A  review  of  59  recent 
overhauls  of  a  variety  of  COMCRUDESLANT  destroyer-type  ships 
indicated  an  approximate  four  month  overhaul  length,  on  the 
average.   The  average  time  between  overhauls  was  approximately 
4  0  months.   For  this  model  therefore,  unless  otherwise  stated, 
the  basic  unit  of  time  is  assumed  to  be  four  months,  although 
the  actual  length  of  the  increment  should  in  general  be  of 
little  importance  to  the  applicability  of  the  model.   Four 
months  should,  however,  provide  sufficient  time  for  the  deci- 
sion maker  to  perceive  a  change  in  condition.   It  provides, 
for  a  ship  life  of,  say,  20  years,  some  60  different  ages. 
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Transition  probabilities  will  be  described  by  p(k|i,t), 
which  is  the  probability  of  the  ship's  being  in  condition  k  at 
age  t+1  given  that  it  was  in  condition  i  at  age  t.   What 
results,  then, is  a  transient  Markov  process  lasting  the  dura- 
tion of  the  ship's  life  (T  time  units)  in  which  a  state  has 
zero  probability  of  occupancy  more  than  once. 

With  this  model,  multi-step  transitional  probabilities 
may  be  calculated.   If  some  initial  state  is  fixed,  the  multi- 
step  transition  probabilities  will  provide  the  probability  of 
entering  a  given  state  (i.e.,  a  given  condition  at  a  given 
age) .   Multi-step  transition  probabilities  will  be  described 
by  P(k,t|k  ,t  )  which  is  interpreted  as  Pr {entering  state 
(k, t) |process  began  in  state  (k  ,t  )}.   Since  at  each  transi- 
tion t  must  increase  by  exactly  1,  the  number  of  steps  in  the 
transition  is  given  simply  by  t-t  .   For  this  model,  multi- 
step  transition  probabilities  are  calculated  by  the  recursion 
formula: 

2 
P(k,t|kQ,to)  =  I   P(i,t-l|k  ,to)p(k|i,t-l), 
i=l 

k  =  1,2. 

Of  particular  interest  in  this  model  will  be  the  probability 
of  entering  a  given  state,  given  that  the  ship  is  new  as  an 
initial  condition.   If  "new"  is  assumed  to  mean  that  the  ship 
is  in  condition  1  at  time  zero  (in  other  words  that  the  ini- 
tial state  is  (1,0)),  then  this  probability  would  be  described 


23 

by  P(k,t|l,0)  for  the  probability  of  entering  state  (k,t). 

For  ease  of  notation  henceforth  this  will  be  written  simply 
as  P (k, t)  . 

This  model  will  provide  multistep  transition  probabili- 
ties (as  has  been  described)  and  therefore  the  probability  of 
the  ship's  being  in  overhaul  at  any  given  age.   From  this  the 
expected  number  of  overhauls  over  the  ship's  life  may  be 
determined.   This  is  simply  the  probability  of  being  in  over- 
haul in  each  time  unit  of  a  ship's  life,  summed  over  the 
ship's  life.   Since  we  are  assuming  that  a  ship  found  to  be  in 

condition  2  will  automatically  undergo  overhaul,  we  may  write: 

T-l 
E (number  of  overhauls  over  ship's  life)  =  J  P(2ft).   The 

t=l 
upper  limit  of  T-l  instead  of  T  reflects  the  assumption  that 

no  ship  would  be  overhauled  during  the  last  time  unit  of  its 

life,  just  in  time  to  be  removed  from  service. 

If  we  desire  to  determine  the  probability  distribution 

of  the  number  of  overhauls  the  ship  will  experience  during 

its  life  (or  up  to  a  given  point  in  its  life) ,  we  must  do  the 

following:   Define  Tr(k,t,x)  as  the  probability  of  the  ship's 

being  in  state  (k,t)  and  having  experienced  x  overhauls 

(including  the  current  one,  if  k=2)  given  that  the  ship  was 

new  at  t=0.   This  probability  may  also  be  calculated  by  means 

of  a  recursion  formula: 
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Tr(k,t,x)  = 


2 

J  ir(i,t-l,x)p(k|i,t-l),k«l, 
i=l 

2 

J  7r(i,t-l,x-l)p(k|i,t-l)  ,k=2. 

i=l 


We  may  then  find  the  probability  distribution  for  the  number 
of  overhauls  in  a  ship's  lifetime  as  follows: 

2 

Pr{ship  will  experience  =  y  (k  T_1      , 

x  overhauls  in  its  lifetime}   ,  ^,    '    r   ' 

k=l 

where  T  is  the  ship's  life  in  terms  of  the  basic  time  unit. 
We  again  use  the  value  for  t=T-l  instead  of  T. 

Now  consider  a  fleet  of  identical  ships,  all  of  the  same 
age.   This  is  obviously  an  oversimplification,  but  will  serve 
to  demonstrate  the  manner  in  which  the  model  may  be  used  in 
considering  groups  of  ships.   Given  the  above  assumptions  and 
independence  between  failures  of  different  ships,  the  proba- 
bility of  having  a  certain  number  of  ships  in  overhaul  at  any 
given  time  out  of  a  fleet  of  a  given  size  follows  a  binomial 
probability  mass  function.   If  a  denotes  the  probability  of 
a  single  ship  being  in  overhaul  at  a  given  age,  the  proba- 
bility of  having  y  ships  of  that  age  in  overhaul  out  of  a 
fleet  size  of  N  ships  is  given  by: 


N 
IYJ 


a-*  (1-a)  J 
This  is  also,  of  course,  the  probability  of  having  N-y  ships 
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available  for  operations  out  of  a  fleet  of  N  ships  (assuming 
that  a  ship  not  in  overhaul  is  available  for  operations) . 

From  this  probability  mass  function  can  be  determined  the 
fleet  size  required  to  provide  some  number  of  operational 
ships  (as  a  minimum)  with  a  specified  probability.   Given  such 
a  requirement,  a  fleet  size  necessary  to  support  it  is  as- 
sumed.  The  probability  of  having  at  least  the  minimum  number 
of  operational  ships  is  computed  using  the  binomial  probabi- 
lity mass  function.   If  this  probability  is  greater  than  the 
specified  probability,  then  the  assumed  fleet  size  is  de- 
creased until  the  requirement  is  just  met.   If  this  proba- 
bility is  less  than  that  specified,  the  assumed  fleet  size  is 
increased  until  the  requirement  is  met. 

Since  the  expected  number  of  overhauls  over  a  ship's  life 
may  be  determined  from  the  model,  assume  that  both  ship  cost 
and  overhaul  cost  are  known.   Then  the  expected  total  cost 
may  be  determined  for  a  single  ship  by  multiplying  the  over- 
haul cost  by  the  expected  number  of  overhauls  and  adding  the 
result  to  the  ship  cost.   Once  this  cost  has  been  calculated 
for  a  single  ship  it  is  known  for  a  fleet  (of  a  given  size)  of 
such  ships. 
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C.   AN  EXAMPLE: 


Before  proceeding  with  development  of  more  comprehensive 
models  it  may  be  helpful  to  demonstrate  application  of  the 
relatively  simple  model  developed  thus  far  in  a  hypothesized 
example. 

Suppose,  for  the  example,  that  the  life  of  the  sample 
ship  is  expected  to  be  10  years,  or  30  time  units;  further, 
suppose  that  the  probability  of  "failure"  increases  in  stages 
as  the  ship  ages,  as  follows: 


p(2|l,t)  = 


.1 
.2 

.3 

.4 
.5 


t=0  , . . 

.,5 

t=6,. . 

.,11 

t=12,. 

..,17 

t=18, . 

.  .  ,23 

t=24,. 

..,29 

It  is  assumed  that  the  ship  begins  in  state  (1,0),  and  that 

if  it  should  "fail",  overhaul  of  the  ship  would  begin 

immediately  and  would,  with  certainty,  return  the  ship  to 

condition  1  at  the  beginning  of  the  next  time  unit. 

The  following  have  been  computed  for  the  example:   Values 

of  P(k,t)  for  all  k  and  t;  values  of  P(k,t|k  ,t  )  for  t 

1  o  o 

through  ten  time  units;  values  of  Tr(k,t,x)  for  all  possible 
k,  t  and  x;  and  the  distribution  and  expected  number  of  over- 
hauls the  ship  will  experience  in  its  life.   All  possible 


P(k,t|k  , t  )  were  not  computed  as  the  result  becomes  somewhat 
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unwieldy.   For  example,  computing  all  values  through  30  time 

units  results  in  approximately  nine  times  the  output  as  for 
ten  units.   Because  of  this,  and  having  demonstrated  the 
manner  in  which  multi-step  transition  probabilities  may  be 
calculated,  henceforth  only  those  values  specifically 
required  for  some  application  will  be  computed,  except  for 
the  special  case  of  P(k,t)  which  is  in  fact  the  family  of 
state  probabilities,  given  that  the  ship  began  in  the  "new" 
condition. 

The  computer  program  used  for  the  calculations  was 
developed  from  the  recursion  relationships  of  the  model.   The 
program,  together  with  the  results  for  the  example,  comprise 
Appendix  A  to  the  thesis.   The  program  is  written  in  the  PL/1 
programming  language  but  the  logic  should  present  no  problem 
to  FORTRAN  programmers . 

The  expected  number  of  overhauls  the  sample  ship  will 
experience  over  its  life  is  6.38.   The  probability  distribu- 
tion of  number  of  overhauls  is  shown  in  Table  1.   These  prob- 
abilities are  determined  by  summing  tt(1,29,x)  and  tt(2,29,x) 
for  all  values  of  x.   A  ship  age  of  29  instead  of  30  is  used 
since  it  is  assumed  that  it  would  be  a  waste  of  resources  to 
overhaul  a  ship  during  the  final  time  unit  of  its  life.   The 
special  family  of  transition  probabilities  P(k,t)  are  given 
by  Table  2.   Since  these  give  the  probability  of  being  in  a 
specific  state,  given  that  the  ship  started  "new",  they  will 
henceforth  be  referred  to  simply  as  "state  probabilities"  to 
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Number  of  Overhauls 
in  Ship's  Life 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 


Probability 

.00002 

.00053 

.00499 

.02287 

.08295 

.17352 

.24346 

.23183 

.14980 

.06494 

.01844 

.00328 

.00034 

.00001 

.00000+ 

.00000+ 


Note: 

Because  of  the  certainty  of  returning  to  condition  1 
after  overhaul,  the  ship  cannot  experience  "failure"  (and 
therefore  overhaul)  during  two  consecutive  time  units.   There- 
fore the  probability  of  more  than  15  overhauls  is  zero. 


TABLE  1 

Probability  Distribution  of  Number  of 

Overhauls  for  the  Two  Condition, 

Age  Dependent  Model  Example 
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P(k,t) 


1 

1 

.900 

2 

1 

.100 

1 

2 

.910 

2 

2 

.090 

1 

3 

.909 

2 

3 

.091 

1 

4 

.909 

2 

4 

.091 

1 

5 

.909 

2 

5 

.091 

1 

6 

.909 

2 

6 

.091 

1 

7 

.818 

2 

7 

.082 

1 

8 

.836 

2 

8 

.164 

1 

9 

.833 

2 

9 

.167 

1 

10 

.833 

2 

10 

.167 

1 

11 

.833 

2 

11 

.167 

1 

12 

.833 

2 

12 

.167 

1 

13 

.750 

2 

13 

.250 

1 

14 

.775 

2 

14 
TABLE  2 

.225 

State  Probabilities  for  the  Two 
Condition,  Age  Dependent 
Model  Example 
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P(k,t) 


1 

15            .767 

2 

15            .233 

1 

16            .770 

2 

16            .230 

1 

17            .769 

2 

17            .231 

1 

18            .769 

2 

18            .231 

1 

19            .692 

2 

19            .308 

1 

20            .723 

2 

20            .277 

1 

21            .711 

2 

21            .289 

1 

22            .716 

2 

22            .284 

1 

23            .714 

2 

23            .286 

1 

24            .715 

2 

24            .285 

1 

25            .643 

2 

25            .357 

1 

26            .679 

2 

26            .321 

1 

27            .661 

2 

27            .339 

1 

28            .670 

2 

28            .330 

1 

29            .665 

2 

29            .335 

1 

30            .667 

2 

30            .333 

TABLE  2  (continued) 
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differentiate  them  from  other  multi-step  transition  proba- 
bil ities . 

Under  the  assumption  that  a  ship  found  to  be  in  condition 
2  will  automatically  undergo  overhaul,  the  state  probability 
P(2,t)  is  also  the  probability  of  the  ship's  being  in  overhaul 
at  age  t.   For  the  present  example,  assume  that  the  decision 
maker  desires  to  maintain  a  fleet  of  ships  identical  to,  and 
of  the  same  age  as,  the  sample  ship.   Assume  further  that  he 
wishes  to  have  a  probability  (say  .95)  of  having  a  certain 
number  of  ships  operational  (that  is,  not  in  overhaul),  say 
30  ships.   The  problem  is  to  determine  the  size  of  fleet 
required  to  give  a  probability  of  at  least  .95  of  30  or  more 
operational  ships . 

With  these  assumptions,  the  required  fleet  size  may  be 
determined  as  follows:   First,  determine  the  highest  probabil- 
ity of  being  in  overhaul  (that  is,  the  maximum  value  of  P(2,2) 
at  any  point  in  the  ship's  life.   For  the  example  this  is 
P(2,25)  (from  Table  2)  which  has  a  value  of  .357.   This  repre- 
sents the  worst  possible  situation.   Next,  using  the  binomial 
probability  mass  function,  obtain  the  probability  distribution 
for  the  number  of  operational  ships  using  the  value  .357  and 
an  assumed  fleet  size.   Increase  or  decrease  the  fleet  size  as 
required  until  the  requirement  of  a  .95  probability  of  30 
operational  ships  is  just  met.   This  is,  of  course,  the  smal- 
lest fleet  size  that  will  meet  the  specified  requirement. 
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The  minimum  fleet  size  satisfactory  to  meet  the  require- 
ments specified  in  the  example  is  56.   The  computer  program 
developed  to  accomplish  this  is  contained  in  Appendix  B.   The 
initial  fleet  size  assumed  was  50  ships.   The  probability  den- 
sity for  the  number  of  operational  ships  for  a  fleet  size  of 
56  is  provided  in  Table  3.   As  would  be  expected,  the  proba- 
bility of  30  or  more  operational  ships  is  .963  (slightly 
greater  than  .95  due  to  discreteness). 

It  should  be  emphasized  that  use  of  the  binomial  proba- 
bility mass  function  is  permitted  only  by  the  fact  that  the 
ships  are  identical  and  of  the  same  age.   The  problem  becomes 
more  complex  (in  computation,  although  not  in  conceptualiza- 
tion) for  a  mixed  fleet. 

Obviously,  with  the  expected  number  of  overhauls  deter- 
mined, if  a  cost  per  overhaul  can  be  estimated  or  is  known, 
then  overhaul  costs  over  a  ship's  life  can  be  calculated. 
Assume  the  following  representative  costs  for  the  example: 
Ship  cost,  $40  million;  overhaul  cost,  $1.5  million.   Then 
the  mean  of  the  total  cost  for  the  example  is: 

$40.0  million  +  6.38  x  $1.5  million 
=  $49.57  million. 

This  can  be  expanded  to  include  the  entire  fleet. 

Although  the  data  is  hypothetical  and  this  model  is  rela- 
tively simple,  the  preceding  example  should  serve  to  indicate 
the  basic  mechanics  and  some  potential  uses  for  such  models. 
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Number  of 

Operational  Ships  Probability 

0  .00000+ 

1  .00000+ 

2  .00000+ 

3  .00000+ 

4  .00000+ 

5  .00000+ 

6  .00000+ 

7  .00000+ 

8  .00000+ 

9  .00000+ 

10  .00000+ 

11  .00000+ 

12  .00000+ 

13  .00000+ 

14  .00000+ 

15  .00000+ 

16  .00000+ 

17  .00000+ 

18  .00000+ 

19  .00000+ 

20  .00001 

21  .00003 

22  .00008 

23  .00022 

24  .00053 

25  .00123 

26  .00264 

TABLE  3 

Probability  Distribution  of  the  Number 
of  Operational  Ships  for  a  Fleet  of 

56  Ships  (Two  Condition,  Age 
Dependent  Example) 
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Number  of 

Operational  Ships  Probability 

27  .00527 

28  .00982 

29  .01706 

30  .02762 

31  .04167 

32  .05857 

33  .07664 

34  .09327 

35  .10548 

36  .11070 

37  .10766 

38  .09685 

39  .08042 

40  .06149 

41  .04317 

42  .02774 

43  .01625 

44  .00864 

45  .00414 

46  .00178 

47  .00068 

48  .00023 

49  .00007 

50  .00002 

51  .00000+ 

52  .00000+ 

53  .00000+ 

54  .00000+ 

55  .00000+ 

56  .00000+ 


TABLE  3  (continued) 
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The  following  chapter  will  deal  with  development  of  more  com- 
prehensive models  based  on  this  one. 


36 
III.   DEVELOPMENT  OF  MORE  COMPREHENSIVE  MODELS 

A.   THE  EFFECT  OF  OVERHAUL  ON  TRANSITION  PROBABILITIES: 

As  a  first  step  in  developing  a  more  comprehensive  model, 
consider  the  following:   If  one  is  willing  to  agree  that  some 
policy  of  overhauling  ships  on  a  regularly  scheduled  basis  may 
be  the  most  economical,  one  is  implicity  assuming  that  an 
overhaul  may  in  some  way  "improve"  succeeding  transition  prob- 
abilities.  Otherwise,  if  no  such  improvement  can  be  attained 
by  an  overhaul,  there  is  no  reason  for  overhauling.   If  a  ship 
is  just  as  likely  to  "fail"  immediately  following  an  overhaul, 
at  some  given  age,  as  it  is  without  a  recent  overhaul,  at  the 
same  age,  then  there  is  obviously  no  justification  for  expen- 
ding the  resources  to  overhaul  the  ship  until  it  actually 
fails . 

It  therefore  appears  that  to  describe  adequately  states 
and  transition  probabilities  requires  knowing  the  time  (i.e., 
number  of  time  units)  since  the  last  overhaul,  as  well  as  the 
ship's  age  and  condition.   Therefore,  let  the  state  indices 
now  be  age  t,  condition  k,  and  time  since  last  overhaul  j; 
the  state  is  described  by  the  ordered  triplet: 

(k,t,j)  =  ship  is  t  units  old  and  has  gone  j  units 
since  its  last  overhaul,  and  is  in  con- 
dition k,  keK,  teH,  jeJ, 
where  again,  K  =  {1,2}  for  a  two  condition  model, 
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H  =  {1,2,...T},  T  being  the  life  of  the  ship 
in  terms  of  the  basic  time  unit, 
and  J  =  {1,2, ...m},  m  being  the  longest  period 
a  ship  could  possibly  go  without  an 
overhaul;  for  a  "demand"  dependent 
(i.e.,  overhaul  only  upon  ship's  reach- 
ing condition  2)  overhaul  policy,  m=T. 

For  a  ship  that  has  yet  to  undergo  its  first  overhaul,  j  will 
be  taken  as  equal  to  t. 

Transition  probabilities  (one-step)  will  be  described  by 
p  (k | i , t, j ) ,  which  is  the  probability  of  the  ship's  being  in 
condition  k  at  age  t+1,  j+1  units  since  its  last  overhaul, 
given  that  it  was  in  condition  i  at  age  t,  j  units  since 
last  overhaul. 

It  will  now  be  necessary  to  describe  multi-step  transi- 
tion probabilities  as  follows: 

p<k,t,j|k0,to,jo), 

Pr{entering  state  (k, t, j ) | transition  began  in 

state  (k  ,  t  ,j  )}. 
o  o  o 

Following  the  convention  adopted  earlier,  the  special  family 
of  transition  probabilites  P  (k, t, j 1 1, 0, 0) ,  the  "state"  proba- 
bilities, will  be  designated  by  P(k,t,j). 

The  expression  7T(k,t,j,x)  will  be  defined  as  the  proba- 
bility of  the  ship's  being  in  state  (k,t,j)  having  experienced 
x  overhauls.   Recursion  formulas  may  be  determined  for  multi- 
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step  transition  probabilities  and  values  of  7r(k,t,j,x),  given 

that  the  process  started  in  state  (1,0,0),  for  this  two-condi- 
tion, age  and  time-since-overhaul  dependent  model.   These 
have  been  developed  in  a  manner  parallelling  those  developed 
earlier  and  are  shown  in  Appendix  C. 
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B.   THE  MULTI -CONDITION,  AGE  DEPENDENT  MODEL: 

While  a  structure  which  provides  for  only  two  categories 
of  material  condition  serves  sufficiently  well  to  demonstrate 
the  mechanics  of  the  model,  to  say  that  a  ship  is  either  in 
good  condition  or  is  in  bad  condition  and  therefore  in  need 
of  overhaul,  with  no  middle  ground,  is  intuitively  unappeal- 
ing.  It  seems  more  likely  that  the  decision  maker  could  group 
his  ships  into  three,  four,  or  even  more  approximately  dis- 
crete groups  based  on  his  subjective  (but  informed)  opinion. 

At  this  point  it  will  be  assumed  that  a  four  condition 
model  is  reasonable  in  that  it  provides  a  wider  range  for 
categorizing  a  ship's  condition  without  providing  so  many 
choices  as  to  make  such  categorization  overly  difficult  and 
arbitrary.   There  will  be  borderline  cases,  to  be  sure,  re- 
gardless of  the  number  of  categories  chosen. 

First,  we  will  consider  a  four  condition  model  in  which 
transition  probabilities  are  unaffected  by  the  time  since 
the  ship  was  last  overhauled  (although  they  still  can  vary 
with  ship  age).   In  this  model,  the  beneficial  effect  of  an 
overhaul  can  be  represented  by  a  change  to  a  better  condition, 
e.g.,  condition  1  with  certainty  or  condition  1  and  condition 
2  with  probabilities  summing  to  unity,  etc.   Transition  prob- 
abilities can  be  used  to  describe  a  gradually  deteriorating 
material  condition;  all  that  is  necessary  is  to  structure  the 
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transition  probabilities  so  that,  over  time,  the  ship's  con- 
dition is  more  likely  to  get  worse  than  it  is  to  get  better. 

With  this  model,  calculation  of  the  entire  set  of  (multi- 
step)  transition  probabilities  remains  reasonable,  although 
still  cumbersome;  further,  we  can  still  determine,  within 
reasonable  computational  limits,  the  probability  distribution 
for  the  number  of  overhauls  which  the  ship  will  experience 
in  its  life,  as  well  as  the  expected  value  of  this  number. 

Computation  of  these  probabilities  is  accomplished  in 
a  manner  identical  to  that  described  in  Chapter  II,  except 
that  the  set  of  possible  conditions  which  the  ship  might 
reach,  K,  is  {1,2,3,4}. 

Such  a  model  will  permit  us  to  make  comparisons  between 
a  demand  dependent  overhaul  policy  and  one  in  which  the  ship 
is  overhauled  at  regularly  scheduled  intervals  (i.e.,  at  cer- 
tain regularly  spaced  points  in  its  life) .   Such  a  scheduled 
overhaul  policy  will  have  the  advantage  of  permitting  at  least 
part  of  the  overhauls  in  a  ship's  life  to  be  scheduled  in 
advance  rather  than  resulting  from  the  ship'  reaching  a  demand 
condition  (i.e.,  unscheduled).   However,  with  this  particular 
model,  we  cannot,  of  course,  base  a  ship's  entry  into  a 
scheduled  overhaul  on  the  time  since  it  was  last  overhauled. 
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C.   THE  MULTI-CONDITION,  AGE  AND  TIME-SINCE-OVERHAUL 
DEPENDENT  MODEL: 

The  most  comprehensive  model  we  will  consider  combines 
both  the  dependence  on  ship  age  and  time  since  overhaul  of 
the  model  of  section  A  with  the  multi-condition  capability  of 
section  B.   This  allows  ship  failure  behavior  to  depend  in  an 
arbitrary  manner  on  ship  age,  time  since  overhaul,  and  present 
material  condition.   This  gain  in  generality  is  not  without 
its  sacrifice.   As  we  will  see,  by  requiring  three  indices 
to  describe  a  state  we  have  lost  the  ability  to  obtain  some 
information  about  the  process  due  to  computational  limita- 
tions.  The  recursion  formulas  for  this  multi-condition,  ship 
age  and  time-since-overhaul  dependent  model  are  shown  in 
Appendix  C. 
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D.   AN  EXAMPLE  FOR  THE  MOST  COMPREHENSIVE  MODEL: 

It  follows  from  the  description  of  this  model  that  there 
are  T  x  m  sets  of  transition  probabilities,  since  probabil- 
ity of  transition  is  dependent  on  both  ship  age  and  time  since 
last  overhaul.   If  one  considers  a  ship's  life,  T,  to  be  about 
25  years,  in  the  case  where  a  ship  isn't  overhauled  until  the 
"demand"  condition  is  reached,  m  is  also  25  years.   This  rep- 
resents 75  time  units  in  the  model,  which  implies  75  x  75  = 
5625  possible  sets  of  transition  probabilities.   This  is 
clearly  unmanageable.   The  answer,  therefore,  is  to  attempt 
to  group  ships  which  seem  to  follow  similar  "failure"  rates 
into  categories  based  upon  age  and  time  since  overhaul. 

Following  this  approach,  use  of  the  multi-condition, 
age  and  time-since-overhaul  dependent  model  will  be  demon- 
strated utilizing  transition  probabilites  estimated  from  the 
RAV  (restricted  availability)  history  of  33  CRUDESLANT  des- 
troyer-type ships.   The  manner  in  which  these  probabilities 
were  estimated  is  described  in  Appendix  D.   The  ships  were 
grouped  into  three  categories  by  age:   Category  I,  1  to  3 
years;  category  II,  4  to  13  years,  and  category  III,  14  to  25 
years,  inclusive.   Three  categories  of  time  since  overhaul 
were  decided  on:   one  (four  month)  time  unit,  two  to  three 
time  units,  and  four  or  more  time  units.   This  results  in  nine 
sets  of  transition  probability  matrices.   These  are  shown  in 
Table  4 . 
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Age  Category  I 
3  =  I 

p(l|l,t,l)  =  0.505 

p(2|l,t,l)  =  0.253 

0)3|l,t,l)  =  0.121 

p(4|l,t,l)  =  0.121 

p(i|k,t,l)  are  of  no  concern  for  i  =  2, 3, 4,  since  ship  can 
only  be  in  condition  1  at  j  =1. 

J  =  2,3 

p(l|l,t,j)  =  0.663  p(l|2,t,j)  =  0.263  p(l|3,t,j)  =  0.263 

p(2|l,t,j)  =  0.263  p(2|2,t,j)  -  0.663  p(2|3,t,j)  =  0.263 

p(3|l,t,j)  =  0.037  p(3|2,t,j)  =  0.037  p(3|3,t,j)  =  0.437 

p(4|l,t,j)  =  0.037  p(4|2,t,j)  =  0.037  p(4|3,t,j)  =  0.037 

j  =  4  or  greater 

p(l|l,t,j)  =  0.682  p(l|2,t,j)  =    0.282  p(l|3,t,j)  =  0.282 

p(2|l,t,j)  =  0.282  p(2|2,t,j)  =    0.682  p(2|3,t,j)  =  0.282 

p(3|l,t,jj  =  0.018  p(3|2,t,j)  =   0.018  p(3|3,t,j)  =  0.418 

p(4|l,t,j)  =  0.018  p(4|2,t,j)  =    0.018  p(4|3,t,j)  =  0.018 


TABLE  4 

One-Step  Transition  Probabilities  for  the 
Multi-Condition,  Age  and  Time-Since- 
Overhaul  Dependent  Example 
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Age  Category  II 

p(l|l,t,l)  =    0.333 

p(2|l,t,l)  =    0.167 

p(3|l,t,l)  =    0.250 

p(4    l,t,l)  -    0.250 

J    =    2,3 

p(l|l,t,j)  =  0.662  p(l|2,t,j)  =  0.262  p(l|3,t,j)  =  0.262 

p(2|l,t,j)  =  0.262  p(2|2,t,j)  =  0.662  p(2|3,t,j)  =  0.262 

p(3|l,t,j)  =  0.038  p(3|2,t,j)  =  0.038  p(3|3,t,j)  =  0.438 

p(4|l,t,j)  =  0.038  p(4|2rt,j)  =  0.038  p(4|3,t,j)  =  0.038 

j  =  4  or  greater 

p(l|l,t,j)  =    0.594  p(l|2,t,j)  =  0.194  p(l|3,t,j)  =  0.194 

p(2|l,t,j)  =    0.194  p(2|2,t,j)  =  0.594  p(2|3,t,j)  =  0.194 

p(3|lrt,j)  =    0,106  p(3|2,t,j)  =  0.106  p(3|3,t,j)  =  0.506 

p(4|lft,j)  =    0.106  p(4|2,t,j)  =  0.106  p(4|3,t,j)  =  0.106 


TABLE  4  (continued) 
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Age   Category   III 

i  =  i 

p(l|l,t,l)  =    0.571 

p(2|l,t,l)  =    0.285 

p(3|l,t,l)  =    0.072 

p(4|l,t,l)  =    0.250 

J    =    2,3 

p(l|l,t,j)  =  0.664  p(l|2,t,j)  =  0.264  p(l|3,t,j)  =  0.264 

p(2|l,t,j)  =  0.264  p(2|2,t,j)  =  0.664  p(2|3,t,j)  =  0.264 

p(3|l,t,j)  =  0.036  p(3|2,t,j)  =  0.036  p(3|3,t,j)  =  0.436 

p(4|l,t,j)  =  0.036  p(4|2,t,j)  =  0.036  p(4|3,t,j)  =  0.036 

j  =  4  or  greater 

p(l|l,t,j)  =    0.627  p(l|2,t,j)  =    0.227  p(l|3,t,j)  =  0.227 

p(2|l,t,j)  =    0.227  p(2|2,t,j)  =    0.627  p(2|3,t,j)  =  0.227 

p(3|l,t,j)  =    0.073  p(3|2,t,j)  =    0.073  p(3|3,t,j)  =  0.473 

p(4|l,t,j)  =    0.073  p(4|2,t,j)  =    0.073  p(4|3,t,j)  =  0.073 


TABLE  4  (continued) 
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The  computer  model  used  for  the  calculations  is  shown  in 

Appendix  E.   This  model  omits  two  features  of  the  model  used 
in  earlier  calculations  (Appendix  A).   These  are:   (1)   The 
computation  of  multi-step  transition  probabilities,  for 
reasons  discussed  earlier,  and  (2)   The  computation  of 
Tr(k,t,j,x),  i.e.,  the  probability  of  being  in  a  given  state 
and  having  experienced  a  given  number  of  overhauls.   The  com- 
putation of  ir(k,t,j,x)  was  omitted  due  to  the  computer  space 
required  to  accomodate  this  four  dimensional  array.   For  the 
demand  dependent  policy,  where  j  may  take  on  any  value  up  to 
t,  this  requires  4  x  75  x  75  x  38  =  755,000  computer  words, 
beyond  the  core  storage  capacity  of  all  but  the  largest  com- 
puters.  Since  determination  of  the  probability  distribution 
of  the  number  of  overhauls  is  the  only  advantage  to  computing 
tt  (k,  t,  j  ,x)  ,  and  since,  as  we  have  seen,  we  can  compute  the 
expected  value  of  this  random  variable  by  other  methods,  under 
the  circumstances  omission  of  this  calculation  seems  a  reason- 
able trade-off,  provided  we  assume  that  decisions  will  be  risk 
neutral  (i.e.,  based  on  expected  cost).   Given  a  computer  with 
sufficient  capacity,  calculation  of  ir(k,t,j,x)  would  be  no 
more  difficult  conceptually  than  it  was  in  the  earlier  model. 

The  state  probabilities  (P(k,t,j))  resulting  from  the 
model  are  calculated  utilizing  the  transition  probabilities  of 
Table  4  for  a  demand  dependent  overhaul  policy.   Because  of 
the  large  number  of  states,  these  probabilites  have  been 
aggregated  over  all  j  for  a  given  pair,  k,t.   This  results  in 
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the  probability  of  the  ship's  being  in  condition  k  at  age  t 
without  consideration  of  the  time  since  last  overhaul.   This 
probability  will  be  denoted  by  <j>(k,t).   Values  of  cj>(k,t)  for 
this  example  are  displayed  in  Appendix  F. 

The  expected  number  of  overhauls  the  ship  will  experience 
over  its  life  is  5.546.   This  value  was  obtained  in  the  same 
manner  as  that  for  the  example  in  Chapter  II:   The  probability 
of  a  ship's  being  in  the  "demand"  condition  (condition  4,  in 
this  case)  at  any  given  ship  age  is  the  probability  of  being 
in  overhaul  at  that  age.   Since  only  one  time  unit  is 
involved,  this  probability  is  also  the  expected  number  of 
overhauls  that  the  ship  will  undergo  at  that  age.   The 
expected  number  of  overhauls  throughout  the  ship's  life  is 
obtained  by  summing  these  values  over  the  ship's  life,  i.e., 
all  possible  ages. 

In  the  following  chapter,  a  comparison  of  the  demand 
dependent  overhaul  policy  reflected  here  will  be  made  with  a 
policy  of  some  fixed  overhaul  cycle  length. 
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E.   RELAXATION  OF  CERTAIN  ASSUMPTIONS: 

Throughout  the  development  of  the  various  models,  we  have 
retained  certain  assumptions:   That,  under  a  demand  dependent 
policy,  a  ship  would  not  be  overhauled  until  it  reached  the 
worst  condition;  that,  when  it  reached  this  condition,  it 
would  be  overhauled  with  probability  1;  and,  that  this  over- 
haul would  return  the  ship  to  condition  1  with  certainty. 
While  we  will  continue  to  make  these  assumptions,  the  approach 
can  accomodate  relaxing  them  with  little  difficulty.   Since 
being  in  the  worst  condition  would  no  longer  be  tantamount 
to  being  in  overhaul,  some  additional  computation  would  be 
required  to  determine  the  probability  of  being  in  overhaul, 
given  the  ship's  age  and  the  length  of  time  since  its  last 
overhaul.   Further,  since  overhaul  would  no  longer  ensure 
the  ship's  returning  to  condition  1,  P(l,t,l)  would  no 
longer  be  equal  to  unity  (for  all  t) ,  and  one  or  more  of 
P(k,t,l)  would  have  values  greater  than  zero  for  k  =  2,3,4. 
Again,  the  change  would  be  more  computational  then  concep- 
tual . 
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IV.   A  COMPARISON  OF  OVERHAUL  POLICIES 

A.   A  SCHEDULED  OVERHAUL  POLICY: 

We  will  now  use  the  most  comprehensive  model  to  consider 
a  policy  by  which  a  ship  is  overhauled  at  some  specified 
interval,  rather  than  waiting  for  some  demand  condition  to  be 
met  before  overhauling.   This  is  the  policy  currently  in 
effect  in  the  Navy.   To  do  this  requires  very  little  modifi- 
cation to  this  model.   For  an  overhaul  cycle  length  (i.e., 
end-of-overhaul  to  end-of -overhaul)  of  j  time  units,  we  modify 
p(k|i,t,j-l)  as  follows: 

p(k|4,t,j-l)  =1,   k  =  1,2,3; 
p(k|i,t,j-l)  =  0,   k  =  1,2,3;   i  =  1,2,3. 

In  this  manner  we  artificially  "force"  a  demand  condition  on 
the  ship.   The  mechanics  of  the  model  do  not  require  any 
change  at  all.   This  modification  does  reflect  the  assumption 
that,  if  the  ship  should  reach  condition  4  prior  to  the  end  of 
the  cycle,  overhaul  would  take  place  ahead  of  schedule  with 
the  next  cycle  beginning  at  that  time. 

We  may  use  the  multi-condition,  age  and  time-since- 
overhaul  dependent  model,  then,  to  compare  policies  of  over- 
hauling at  some  fixed  interval  with  the  demand  dependent  pol- 
ciy  described  earlier.   We  will  utilize  our  previous  example 
(Section  D  of  Chapter  III)  to  demonstrate  this. 
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With  the  model  modified  to  reflect  a  fixed  overhaul  pol- 
icy as  described  above,  the  results  shown  in  Table  5  were 
obtained.   In  the  example  of  Chapter  II  we  used  the  highest 
probability  of  being  in  overhaul  over  the  ship's  life  as  a 
basis  for  determining  the  necessary  fleet  size  to  meet  speci- 
fied requirements.   This  highest  probability  is  shown  in  Table 
5,  along  with  the  corresponding  value  of  j,  as  well  as  the 
cycle  length  (in  years)  and  the  expected  number  of  overhauls 
the  ship  will  experience  under  that  policy.   The  symbol  "°°" 
implies  a  strictly  demand  dependent  policy.   We  observe  that 
the  highest  probability  of  being  in  overhaul  occurs  at  the 
end  of  the  first  cycle  in  each  case  (i.e.,  at  j  =  3  for  a  one 
year  cycle  length,  j  =  6  for  a  two  year  cycle  length,  etc.). 
To  use  these  values  to  determine  a  required  fleet  size  would 
imply  a  whole  fleet  of  exactly  the  same  age.   Although  we 
accepted  this  assumption  in  the  example  of  Chapter  II,  such 
an  assumption  is  obviously  unrealistic  and  would  distort  the 
decision  making  process  in  this  case.   Even  if  all  the  ships 
in  a  fleet  were  new  at  the  same  time,  we  must  realize  that  in 
actuality  the  decision  maker  would  stagger  the  point  at  which 
ships'  overhaul  cycles  were  begun.   Consequently,  we  must  find 
some  way  to  overcome  this  problem. 
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Overhau 

Cycle  Len 

(Years) 

1 
gth 

Expected 
Number  of 
Overhauls 
(Lifetime) 

Highest 
Probability 
of  Being  in 

Overhaul 

Ship  Age  at 
which  Highest 
Probability 
Occurs 

1 

25.67 

.879 

J  =   3 

2 

13.96 

.811 

j  =   6 

3 

10.24 

.765 

j  =   9 

4 

8.43 

.641 

j  =  12 

5 

7.42 

.480 

j  =  15 

6 

6.82 

.371 

j  =  18 

7 

6.42 

.294 

j  =  21 

8 

6.16 

.238 

j  =  24 

9 

5.98 

.199 

j  =  27 

10 

5.85 

.171 

j  =  30 

5.55 


.121 


j  =   2 


where  j  is  the  ship's  age  in  terms  of  the  basic  time  unit. 


TABLE  5 


Expected  Number  of  Overhauls  and  Highest 
Probabilities  of  Being  in  Overhaul  for 
Various  Overhaul  Policies 
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B.   THE  MIXED  FLEET: 

It  is  much  more  likely  that  in  reality  only  a  few  new 
ships  would  be  acquired  at  a  time.   We  will  assume  that  our 
ship  acquisition  policy  calls  for  purchasing  some  number  of 
new  ships  every  two  years.   If  we  assume  that  this  takes  place 
over,  say,  ten  years,  then  we  essentially  have  five  different 
types  of  ship  to  consider,  since  the  probability  of  being  in 
overhaul  varies  with  ship  age,  even  though  we  still  assume 
that  the  ships  are  otherwise  identical. 

In  actual  practice  in  the  Navy  a  fleet  of  ships  is  con- 
tinuously experiencing  replacement  of  old  ships  with  new  ones. 
Consequently  we  should  examine  our  "mixed"  fleet  over  some 
period  where  we  have  essentially  reached  a  "steady  state"  in 
number  of  ships  (rather  than  some  point  where  we  are  either 
still  acquiring  ships  on  net  or  have  begun  decreasing  the 
total  fleet  size) ,  say  a  five  year  period  beginning  when  our 
ships1  ages  vary  from  12  to  2  0  years.   We  will  wish  to  examine 
the  mixed  fleet  at  perhaps  five  points  in  time  during  this 
period,  say  a  year  apart,  in  order  again  to  determine  the 
total  fleet  size  required  to  provide  us  with  a  given  number  of 
operational  ships  with  a  specified  probability.   Unfortu- 
nately, since  all  our  ships  do  not  simultaneously  have  the 
same  probability  of  being  in  overhaul  we  cannot  apply  the 
binomial  probability  mass  function  as  we  did  in  Chapter  II. 
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There  is,  however,  another  method  of  approaching  the 
problem.   First,  assume  that  we  have  a  relatively  large  num- 
ber of  ships.   Secondly,  suppose  that  the  probability  of  one 
ship's  being  in  overhaul  is  independent  of  the  probability  of 
any  other  ship's  being  in  overhaul.   This  should  in  general 
be  true  except  for  combatant  situations,  natural  disasters, 
etc.,  or  unless  the  capacity  of  the  overhauling  shipyards 
becomes  a  constraint. 

If  these  conditions  are  satisfied  then  we  may  apply  the 
central  limit  theorem  of  probability.   Essentially,  the  cen- 
tral limit  theorem  states  (for  our  application)  that  the  sum 
of  a  sequence  of  N  independently  distributed  random  variables 
converges  in  distribution  as  N-*-°°  to  a  random  variable  that  is 
normally  distributed  with  mean  equal  to  the  sum  of  the  com- 
ponent means  and  with  variance  equal  to  the  sum  of  the  com- 
ponent variances,  provided  that  Liapunov's  condition  is  met. 
Cramdr'*  states  Liapunov's  condition  as  follows:   Let  z,  , 
z2/««-  be  independent  random  variables,  and  denote  by  m.  and 
a^   the  mean  and  the  standard  deviation  of  z. .   Suppose  that 

the  third  absolute  moment  of  z.  about  its  mean 

1 

pi?  =  E(|zi  -  m±\3) 
is  finite  for  every  i,  and  write 


k See  Harald  Cramer,  Mathematical  Methods  of  Statistics, 
Princeton  University  Press,  Princeton,  NJ,  1946,  pp.  215-218. 
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Pl  +  p2  +  P3  +  **•  +  p] 


If  the  condition 


lim     £  =  o 
k  -h»  a 

k 
is  satisfied,  then  the  sum  z  =  £  z .  is  asymptotically  normal 

i=l  1 
(m,a)  where  m  and  a   are  given  by: 


m  =  m,  +  m~  +  . 
a2  =  a2  +  a2  +  • 


.  + 


V 


.  +  a,  . 


For  a  set  of  Bernoulli  trials 


p3  =  E(|zi  -  pj3)  =  p-q^p?  +  q?) 


where 


Therefore 


Now: 


Therefore 


<  p.q. 


q±  -  i  -  p±. 


p3  <  I  Pi<i  -  Pi). 

i=l        x 


a2  =  I  p±(l  -  P±) 

i=l  x      x 


p  < 


i=l 


Pia 


pi> 


-i/« 
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CO 

If  the  series   £  p. (1  -  p.)  is  divergent ,  Liapunov's  condition 

i=l  1      1 
is  satisfied  and  thus  the  variable  z  is  asymptotically  normal 


jP.,/[p.(i  -p.) 

i=l    i=l 

The  series   J  p. (1  -  p.)/  where  p.  is  the  probability  of  the 

i=l  x      X  X 

ith  ship's  being  in  overhaul,  will  be  divergent  if  for  all  i 

p.  equals  neither  0  nor  1.   This  can  be  reasonably  assumed  for 

our  application. 

The  approximating  normal  distribution  will  then  have  a 
N  N 

mean  of  \   p.  and  a  variance  of   Y  p.  (1  -  p.  )  . 
i=l  1  i=l  1      1 

Since  we  have  our  ships  in  five  groups  and  are  assuming 

that  all  ships  within  a  group  are  identical,  we  may  simplify 

the  summation  process  somewhat  as  follows:   If  we  have  n. 

ships  with  probability  p.  of  being  in  overhaul  at  some  speci- 

1  5 

fied  time,  i  =  1,...,5,  then  we  may  write:   y  =  J,   n.p.,  and 

5  i=l  1  1 

s2  =  ][  n.p.  (1  -  p.),  in  order  to  obtain  the  mean,  y,  and  the 

i=l  x   x  1 

variance,  s2,  of  our  approximating  normal  distribution. 

5 
If  we  take  N  =  J  n.,  then,  using  the  normal  distribution, 
i=l  1 
we  may  approximate  the  probability  of  having  k  out  of  N  ships 

in  overhaul,  where  k  is  equal  to  or  less  than  a  specified 

value  b,  as  follows: 


Pr{k  <  b}  :  »[b  +  1/2  '  v] 


$(yQ)  is  defined  as  the  cumulative  distribution  function  for 
the  unit  normal  probability  density  function,  and 
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»(y  )  .  JL  /Y°  e^2/2dy. 
°     /2?  -« 

The  +1/2  term  is  added  to  the  argument  of  $  in  order  to 
reflect  the  discreteness  of  the  distribution  we  are  approxi- 
mating.  The  reasoning  here  is  similar  to  that  in  the 
DeMoivre-LaPlace  limit  theorem. 

Using  this  approach,  then,  we  may  determine  the  size  of 
fleet  necessary  to  meet  our  requirement.   Obviously,  given 
that  we  have  initially  assumed  a  fleet  size  as  before,  any 
increase  or  decrease  may  be  implemented  among  our  five  dif- 
ferent groups  of  ships  in  various  ways.   Where  we  decide  to 
make  such  an  increase  or  decrease  may  have  some  effect  on  our 
outcome . 

Now  that  we  have  developed  the  approach  to  the  problem 
using  the  normal  approximation,  there  is  no  reason  why  we  must 
restrict  ourselves  to  ships  which  are  identical  except  for 
age.   We  are  equipped  to  consider  ships  with  entirely  differ- 
ent sets  of  material  condition  transition  probabilities.   For 
each  type  we  must,  of  course,  go  through  the  procedure  demon- 
strated in  the  example  of  Chapter  III  for  determining  the 
state  probabilities. 
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C.   AN  EXAMPLE: 

Returning  to  the  results  of  the  example  of  Chapter  III: 
Suppose  we  wish  to  consider  a  fleet  of  five  different  groups 
of  such  ships,  initially  varying  in  age  from  twelve  to  twenty 
years,  over  a  five  year  period.   Define  t  as  the  epoch  in  time 
under  consideration.   For  convenience  let  x  be  equal  to  the 
age  of  the  oldest  group  of  ships.   Table  6  shows  values  for  t 
and  the  ages  of  the  different  groups,  in  terms  of  the  basic 
(four  month)  time  unit,  at  each  of  our  five  observations. 

For  the  example  we  will  assume  that  the  decision  maker 
requires  a  probability  of  .95  for  the  event  {20  or  more  ships 
not  in  overhaul}.   We  will  initially  assume  a  required  fleet 
size  of  25  ships:   five  in  each  group.   If  it  is  necessary  to 
increase  the  fleet  size  to  meet  the  requirement,  we  will  begin 
by  adding  a  ship  to  Group  5,  then  one  to  Group  4,  etc.   If,  on 
the  other  hand,  we  wish  to  decrease  the  fleet  size  because  the 
assumed  size  is  more  than  sufficient  to  meet  the  requirement, 
we  will  remove  one  ship  from  Group  1,  then  one  from  Group  2, 
etc. 

In  the  example  we  will  consider  overhaul  cycle  lengths  of 
four,  six,  eight,  and  ten  years  as  well  as  the  demand  depen- 
dent case.   For  each  of  these  we  will  take  the  observation 
with  the  largest  fleet  size  necessary  to  support  the  decision 
maker's  requirement  and  compare  the  policies  on  a  mean  value 
basis,  economizing  total  cost.   Here  we  will  assume,  as  we 
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Age  of 

Group  1  Age  of  Age  of  Age  of  Age  of 

Observation   t    (oldest)   Group  2  Group  3   Group  4  Group  5 

1  60      60  54  48  42  36 

2  63      63  57  51  45  39 

3  66      66  60  54  48  42 

4  69      69  63  57  51  45 

5  72      72  66  60  54  48 

Age  is  given  in  terms  of  the  basic  time  unit. 


TABLE  6 

Ship  Ages  for  the  Mixed 
Fleet  Example 
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have  throughout  the  thesis,  that  our  decision  maker  is  not  a 

risk  averter  and  will  therefore  base  his  decision  on  expected 

cost. 

Table  7  shows  the  values  of  p.  (the  probability  of  being 
in  overhaul  for  the  ith  group)  for  each  of  the  five  observa- 
tions under  each  of  the  five  overhaul  policies.   Using  these 
values  we  now  compute  the  number  of  ships  in  each  group  (under 
the  decision  rules  previously  established)  needed  to  meet  our 
requirement.   The  computer  model  used  to  do  this  is  shown  as 
Appendix  G.   This  model,  in  the  form  shown,  actually  computes 
the  probabilities  for  the  event  {20  or  more  ships  not  in  over- 
haul} for  all  combinations  of  four,  five,  and  six  ships  in 
each  of  our  five  groups,  so  the  decision  maker  could  examine 
various  rules  for  increasing  or  decreasing  the  initially 
assumed  fleet  size. 

The  requirements  resulting  from  the  computation  are  tabu- 
lated in  Table  8  for  each  of  the  five  overhaul  policies  we  have 
chosen  to  consider.   The  n.  are  the  requirement  for  the  number 
of  ships  in  each  of  the  five  groups  with  N  being  total  fleet 
size. 

For  the  various  policies,  the  maximum  N  (upon  which  we 
will  base  our  decision)  is  summarized  as  follows: 
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Four  Year  Overhaul  Cycle 


Observation 

fl 

h. 

El 

!i 

£5 

1 

.11683 

.10314 

.12567 

.09749 

.17841 

2 

.10860 

.10721 

.11084 

.10418 

.11817 

3 

.10600 

.11683 

.10314 

.12567 

.09749 

4 

.10830 

.10860 

.10721 

.11084 

.10418 

5 

.11260 

.10600 

.11683 

.10314 

.12567 

Six  Year  Overhaul  Cycle 


Observation 

^1 

£i 

P_3 

!i 

!s. 

1 

.08313 

.09742 

.08469 

.07609 

.15261 

2 

.08511 

.09094 

.08631 

.08047 

.10574 

3 

.08619 

.08313 

.09742 

.08469 

.07609 

4 

.08656 

.08511 

.09094 

.08631 

.08047 

5 

.08980 

.08619 

.08313 

.09742 

.08469 

TABLE  7 

Probabilities  of  Being  in  Overhaul  at  Each 
Observation  for  the  Mixed  Fleet  Example 
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Eight  Year  Overhaul  Cycle 


Observation 

ll 

!l2 

!i 

!!i 

^5 

1 

.07856 

.07533 

.08715 

.06887 

.10810 

2 

.08183 

.07668 

.07215 

.07217 

.10670 

3 

.07494 

.07856 

.07533 

.08715 

.06887 

4 

.07598 

.08183 

.07668 

.07215 

.07217 

5 

.07875 

.07494 

.07856 

.07533 

.08715 

Ten  Year  Overhaul  Cycle 


Observation 

Pi 

!i 

h. 

!i- 

P-± 

1 

.07678 

.06968 

.06854 

.06591 

.10064 

2 

.07114 

.07034 

.06909 

.06818 

.10061 

3 

.07252 

.07678 

.06968 

.06854 

.06591 

4 

.07466 

.07114 

.07034 

.06909 

.06818 

5 

.07064 

.07252 

.07678 

.06968 

.06854 

TABLE  7  (continued) 
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Demand  Dependent  Policy 


Observation 

£l 

!_2 

^3 

l± 

^1 

1 

.06395 

.06395 

.06394 

.06165 

.09980 

2 

.06395 

.06395 

.06395 

.06407 

.09980 

3 

.06395 

.06395 

.06395 

.06394 

.06165 

4 

.06395 

.06395 

.06395 

.06395 

.06407 

5 

.06395 

.06395 

.06395 

.06395 

.06394 

TABLE  7  (continued) 
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Four  Year  Overhaul  Cycle 


Observation 

nl 

n2 

n3 

n4 

n5 

N 

1 

5 

5 

5 

5 

6 

26 

2 

5 

5 

5 

5 

5 

25 

3 

5 

5 

5 

5 

5 

25 

4 

5 

5 

5 

5 

5 

25 

5 

5 

5 

5 

5 

5 

25 

Six  Year  Overhaul  Cycle 


Observation 

nl 

n2 

n3 

n4 

n5 

N 

1 

5 

5 

5 

5 

5 

25 

2 

4 

5 

5 

5 

5 

24 

3 

4 

5 

5 

5 

5 

24 

4 

4 

5 

5 

5 

5 

24 

5 

4 

5 

5 

5 

5 

24 

TABLE  8 

Fleet  Size  Requirements  for  the 
Mixed  Fleet  Example 
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Eight  Year  Overhaul  Cycle 


Observation 

nl 

n2 

n3 

n4 

n5 

N 

1 

4 

5 

5 

5 

5 

24 

2 

4 

5 

5 

5 

5 

24 

3 

4 

5 

5 

5 

5 

24 

4 

4 

5 

5 

5 

5 

24 

5 

4 

5 

5 

5 

5 

24 

Ten  Year  Overhaul  Cycle 


Observat 

ion 

nl 

n2 

n3 

n4 

n5 

N 

1 

4 

5 

5 

5 

5 

24 

2 

4 

5 

5 

5 

5 

24 

3 

4 

5 

5 

5 

5 

24 

4 

4 

5 

5 

5 

5 

24 

5 

4 

5 

5 

5 

5 

24 

TABLE  8  (continued) 
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Demand  Dependent  Policy 


Observation 

nl 

n2 

n3 

n4 

n5 

N 

1 

4 

5 

5 

5 

5 

24 

2 

4 

5 

5 

5 

5 

24 

3 

4 

4 

5 

5 

5 

23 

4 

4 

4 

5 

5 

5 

23 

5 

4 

4 

5 

5 

5 

23 

TABLE  8  (continued) 
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POLICY  Nmax 

Four  Year  Overhaul  Cycle  26 

Six  Year  Overhaul  Cycle  25 

Eight  Year  Overhaul  Cycle  24 

Ten  Year  Overhaul  Cycle  24 

Demand  Dependent  Policy  24 

Recalling  the  expected  number  of  overhauls  from  Table  5, 
we  may  compute  total  expected  cost  for  the  fleet,  under  each 
policy,  as  follows: 

Total  expected  cost  =  N   [Ship  cost  + 
Overhaul  cost  (Expected  number  of  overhauls)] 

As  before,  we  will  assume  a  ship  cost  of  $40  million  and  an 
overhaul  cost  of  $1.5  million.   Thus  for  the  four  year  overhaul 
cycle  policy,  total  expected  cost  would  be  computed  as  follows: 

Total  expected  cost  =  26[$40  x  106+  $1.5  x  106(8.43)] 

=  $1.37  billion. 

For  the  remaining  policies  considered,  total  expected  cost  is: 

POLICY  TOTAL  EXPECTED  COST 

Six  Year  Overhaul  Cycle  $1.26  billion 

Eight  Year  Overhaul  Cycle  $1.18  billion 

Ten  Year  Overhaul  Cycle  $1.17  billion 

Demand  Dependent  Policy  $1.16  billion 
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The  decision  maker  would  then  choose  the  demand  dependent  pol- 
icy as  that  which  would  be  most  economical  (based  on  expected 
costs).   As  one  would  expect,  and  as  we  see,  due  to  the  fact 
that  ship  cost  is  relatively  large  compared  to  overhaul  cost, 
the  total  expected  cost  is  more  sensitive  to  required  fleet 
size  than  to  expected  number  of  overhauls. 

Suppose  now  that  the  decision  maker  is  willing  to  relax 
his  requirement  to  a  probability  of  .90  for  the  event  {20  or 
more  ships  not  in  overhaul}.   In  the  same  manner  as  before,  we 
obtain  the  following  results: 

POLICY  N     TOTAL  EXPECTED  COST 

Four  Year  Overhaul  Cycle  25  $1.32  billion 

Six  Year  Overhaul  Cycle  24  $1.21  billion 

Eight  Year  Overhaul  Cycle  24  $1.18  billion 

Ten  Year  Overhaul  Cycle  23  $1.12  billion 

Demand  Dependent  Policy  23  $1.11  billion 

On  the  other  hand,  if  he  requires  a  probability  of,  say,  .98 
for  the  event  {20  or  more  ships  not  in  overhaul},  we  obtain 
the  following: 

POLICY  N     TOTAL  EXPECTED  COST 

Four  Year  Overhaul  Cycle  27  $1.42  billion 
Six  Year  Overhaul  Cycle  26  $1.31  billion 
Eight  Year  Overhaul  Cycle  25  $1.23  billion 
Ten  Year  Overhaul  Cycle  24  $1.17  billion 
Demand  Dependent  Policy     24       $1.16  billion 
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We  observe  that,  while  the  decision  is  relatively  insensitive 
to  such  a  change  in  requirements,  the  required  size  of  fleet 
is  not. 

Up  to  this  point  we  have  assumed  that  the  cost  of  a 
scheduled  overhaul  (under  a  policy  calling  for  an  overhaul 
cycle  of  some  specified  length)  is  the  same  as  that  of  an 
unscheduled  overhaul  (i.e.,  one  resulting  from  the  ship's 
reaching  the  demand  condition  prior  to  the  end  of  the  overhaul 
cycle).   For  a  strictly  demand  dependent  policy,  of  course, 
all  overhauls  are  unscheduled.   An  unscheduled  overhaul  should 
be  expected  to  cost  more  than  a  scheduled  one,  due  to  uncer- 
tainty as  to  the  number  of  ships  in  overhaul  at  any  given 
time,  the  inability  to  plan  specific  work  ahead  of  time  based 
on  an  individual  ship,  etc. 

Using  our  multi-condition,  age  and  time-since-overhaul 
dependent  model,  the  expected  number  of  both  scheduled  and 
unscheduled  overhauls  which  a  ship  will  experience,  under  a 
policy  calling  for  specified  cycle  lengths,  can  be  computed  as 
follows: 

T-l 
E (number  of  scheduled  overhauls)  =  £  P(4,t,L  ); 

t=l       C 

T-l  V1 
E  (number  of  unscheduled  overhauls  -   ][     ][  P(4,t,j) 

t=l   j=l 

where  Lc  is  the  specified  overhaul  cycle  length.   These  values 
have  been  computed  for  our  example  and  are  shown  in  Table  9 
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Overhaul 


Cycle 
Length 
(Years) 

Expected  Number 
of  Scheduled 
Overhauls 

Expected  Number 
of  Unscheduled 
Overhauls 

Total 

1 

21.70 

3.97 

25.67 

2 

-9.95 

4.01 

13.96 

3 

5.59 

4.65 

10.24 

4 

3.48 

4.95 

8.43 

5 

2.27 

5.15 

7.42 

6 

1.54 

5.28 

6.82 

7 

1.06 

5.36 

6.42 

8 

0.74 

5.42 

6.16 

9 

0.52 

5.46 

5.98 

10 

0.37 

5.48 

5.85 

0.00 


5.55 


5.55 


TABLE  9 


Expected  Number  of  Scheduled  and 
Unscheduled  Overhauls  for  the 
Multi-Condition,  Age  and  Time- 
Since-Overhaul  Dependent 
Model  Example 
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for  overhaul  cycle  lengths  of  one  to  ten  years  and  for  the 
demand  dependent  policy. 

Suppose  we  assume  that  the  cost  per  overhaul  of  an 
unscheduled  overhaul  is  $2  million,  while  that  of  a  scheduled 
overhaul  is  $1.5  million.   We  will  again  examine  the  total 
expected  fleet  cost,  taking  this  assumption  into  account.   If 
we  require  a  probability  of  .90  for  the  event  {20  or  more 
ships  not  in  overhaul}  (denote  this  as  p    (20)),  we  obtain 
the  following  results: 

POLICY  N 

Four  Year  Overhaul  Cycle  25 

Six  Year  Overhaul  Cycle  24 

Eight  Year  Overhaul  Cycle  24 

Ten  Year  Overhaul  Cycle  2  3 

Demand  Dependent  Policy  23 


TOTAL  EXPECTED  COST 

$1.38  billion 
$1.27  billion 
$1.25  billion 
$1.18+  billion 
$1.18-  billion 


For  a  p    (20)  of  .95  we  obtain: 
rreq 

POLICY  N 

Four  Year  Overhaul  Cycle  26 

Six  Year  Overhaul  Cycle  25 

Eight  Year  Overhaul  Cycle  24 

Ten  Year  Overhaul  Cycle  24 

Demand  Dependent  Policy  24 


TOTAL  EXPECTED  COST 

$1.43  billion 

$1.32  billion 

$1.25  billion 

$1.24  billion 

$1.23  billion 


And,  for  p    (20)  =  .98: 

req 
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POLICY  N     TOTAL  EXPECTED  COST 

Four  Year  Overhaul  Cycle  27  $1.49  billion 

Six  Year  Overhaul  Cycle  26  $1.37  billion 

Eight  Year  Overhaul  Cycle  25  $1.30  billion 

Ten  Year  Overhaul  Cycle  24  $1.24  billion 

Demand  Dependent  Policy  24  $1.23  billion 

As  we  can  see,  the  decision  as  to  the  most  economical  policy 
would  remain  unchanged  with  the  assumption  that  an  unscheduled 
overhaul  costs  $2.0  million  ($0.5  million  more  than  a  sched- 
uled one)  . 

We  may  approach  the  trade-off  between  scheduled  and 
unscheduled  overhauls  from  another  viewpoint,  however.   Sup- 
pose we  wish  to  find  some  coefficient  c,  where 

cost  of  unscheduled  overhaul 
c  = 


cost  of  scheduled  overhaul  ' 


such  that  the  decision  maker  would  be  indifferent  between  two 
alternative  overhaul  policies.   To  illustrate  the  manner  in 
which  we  can  do  this,  consider  the  ten  year  overhaul  cycle 
policy  and  the  demand  dependent  policy  for  p    (20)  =  .95. 
Since  the  required  fleet  size  is  the  same,  we  need  only  con- 
sider single  ship  expected  overhaul  cost  in  order  to  determine 
a  value  for  c.   We  do  this  as  follows: 

E (number  of  scheduled  overhauls  for  10  yr .  policy) 

+  c  x  E (number  of  unscheduled  overhauls  for  10  yr.  policy) 
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=  E (number  of  scheduled  overhauls  for  demand  dependent 
policy)  +  c  x  E (number  of  unscheduled  overhauls  for 
demand  dependent  policy) . 

Substituting  the  values  from  Table  9  we  get: 

.37  +  5.48c  =  .00  +  5.55c 

c  =  ^1   =  5  3 

C      Q7  D.J 

This  implies  that  for  the  particular  set  of  transition  proba- 
bilities assumed  for  the  example,  the  cost  of  an  unscheduled 
overhaul  would  have  to  be  of  the  order  of  five  times  that  of 
a  scheduled  one  in  order  for  the  decision  maker  to  change  his 
decision  to  adopt  the  demand  dependent  overhaul  policy. 

In  our  example  the  demand  dependent  overhaul  policy  has 
consistently  been  the  choice,  economically.   This  will  always 
be  true  unless  we  impute  a  penalty  to  the  cost  of  an  unsched- 
uled overhaul.   Given  such  a  penalty,  and  different  behavior 
of  transition  probabilities  from  our  example  (especially  if 
probability  of  "failure"  continues  to  increase  as  ship  age 
and  time  since  overhaul  increase) ,  we  could  easily  have 
reached  a  different  decision. 

We  have  shown,  then,  how  our  model  can  be  used  to  compare 
alternative  overhaul  policies  by  considering  total  expected 
cost  of  a  fleet  of  ships  in  addition  to  the  cost  differential 
between  scheduled  and  unscheduled  overhauls. 
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V.   RECOMMENDATIONS,  SUMMARY,  AND  CONCLUSIONS 

A.   RECOMMENDATIONS  FOR  FURTHER  WORK: 

The  next  step,  and  the  most  important  toward  actual 
implementation  of  the  models,  and  in  particular  the  most  com- 
prehensive (multi-condition,  age  and  time-since-overhaul 
dependent)  model,  is  determination  of  reasonable  estimators 
for  transition  probabilities  based  on  the  best  data  available. 
Once  these  have  been  determined,  the  model  implemented,  and 
transition  probabilities  refined  so  that  the  results  of  the 
model  are  relatively  reliable,  the  model  may  be  used  to  exa- 
mine a  large  number  of  possible  maintenance  strategies.   For 
example,  the  decision  maker  might  want  to  consider  the  alter- 
native of  assigning  his  ships  to  shipyards  upon  their  reaching 
condition  3  for  limited  maintenance  not  actually  amounting  to 
overhaul  (as  we  have  used  the  term  here) .   There  are  no  doubt 
other  maintenance  policies  of  this  nature  that  can  be  eval- 
uated using  this  multi-condition,  time-since-overhaul  depen- 
dent model  as  a  basis. 

If  it  were  found  to  be  reasonable  to  describe  transition 
probabilities  analytically  in  terms  of  ship  age  and  time  since 
last  overhaul,  the  model  certainly  could  be  modified  to  accom- 
modate such  a  change.   As  a  matter  of  fact,  transition  proba- 
bilities could  be  computed  internally  by  the  computer  model. 
This  would  eliminate  the  necessity  of  reading  in  (in  some 
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manner)  long  lists  of  transition  matrices,  as  we  have  done 
here. 

Efforts  to  develop  criteria  for  an  acceptable  material 
condition  index  for  ships  should  continue,  although  the  com- 
plexity and  variety  of  configurations  encountered  in  warships 
does  not  make  this  an  easy  assignment.   An  index  of  some  kind 
is  essential  if  one  is  to  expect  valid  results  from  this 
model.   In  the  interim,  however,  the  model  should  be  imple- 
mented utilizing  whatever  rough  index  of  material  condition 
that  can  be  developed  in  the  short  run. 

It  must  be  emphasized  that  results  obtained  from  the 
models  are  only  as  good  as  their  input.   The  models  are  geared 
to  categorization  of  material  condition  either  by  subjective 
judgment  or  by  some  quantitative  scheme  for  determining  a 
material  condition  index.   Critical  to  the  model's  usefulness 
as  an  aid  to  decision  making  is  accurate  determination  of 
transition  probabilities,  as  indicated  above.   As  pointed  out 
in  Appendix  D,  these  should  be  estimated  initially  from  the 
best  historical  data  available,  using  criteria  which,  once 
decided  upon,  should  be  applied  consistently  to  all  ships. 
Then,  using  the  Bayesian  approach,5  transition  probabilities 
should  be  updated  as  actual  ship  behavior  is  observed. 


5 For  one  example  of  how  this  approach  can  be  used,  see 
Chapter  5  of  John  W.  Devanney  III,  Marine  Decisions  Under 
Uncertainty,  National  Sea  Grant  Program,  Sea  Grant  Project 
GH-88 ,  M. I ,T. ,  Cornell  Maritime  Press,  Cambridge,  Maryland, 
November  1971. 
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We  have  demonstrated  the  model  using  only  one  set  of 
transition  probability  matrices.   It  would  be  useful  to  hypo- 
thesize transition  probabilities  based  on  different  assumed 
forms  of  behavior  of  the  ship's  material  condition  (e.g., 
deterioration  more  or  less  linear  over  time,  etc.)  and  compare 
the  results  with  those  obtained  here.   In  this  manner  some 
"feel"  may  be  obtained  for  the  dynamics  of  the  model. 

An  approach  to  modeling  the  change  in  a  ship's  condition 
over  time  that  deserves  some  consideration  involves  the  use  of 
a  semi -Markov  process  model.   In  a  semi -Markov  process,  the 
states  occupied  on  successive  transitions  are  governed  by  the 
transition  probabilities  (as  in  our  model)  of  a  Markov  process 
(called  the  imbedded  Markov  process) .   The  time  the  process 
stays  in  a  particular  state,  however,  is  a  random  variable 
(either  discrete  or  continuous)  described  by  an  exponential 
density  (Poisson  process) .   This  random  variable  is  referred 
to  as  the  "holding"  time  and  is  dependent  on  the  state  pres- 
ently occupied  and  on  the  state  to  which  the  next  transition 
will  be  made.   Howard,  Volume  II,6  provides  a  rather  thorough 
discussion  of  semi-Markov  processes. 


6 Ronald  A.  Howard,  Dynamic  Probabilistic  Systems,  Volume 
II :  Semi-Markov  and  Decision  Processes ,  John  Wiley  and  Sons, 
New  York,  1971,  Chapters  10  and  11. 
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B.   SUMMARY  AND  CONCLUSIONS: 

Beginning  with  a  very  simple  description  of  the  behavior 
of  a  ship's  material  condition,  we  have  developed  a  set  of 
models  which,  while  still  conceptually  simple,  is  capable  of 
handling  a  rather  wide  range  of  ship  deterioration  behavior. 
The  approach  used  in  developing  the  models  considers  the  ship 
as  a  single  component,  rather  than  a  set  of  many  smaller  com- 
ponents.  Modeling  the  behavior  of  the  ship's  material  con- 
dition over  time  by  a  transient  Markov  process  represents  an 
unusual  if  not  unique  approach  to  this  area  of  no  small  con- 
cern to  the  Navy. 

With  the  central  limit  theorem,  as  we  have  seen,  there  is 
no  problem  in  going  from  considerations  of  a  fleet  of  identi- 
cal ships,  all  of  the  same  age,  to  a  "mixed"  fleet,  where  not 
only  are  ships  not  all  of  the  same  age  but  are  dissimilar  in 
other  respects  as  well. 

We  have  seen  how  the  models  can  be  used  to  determine  the 
effect  of  various  overhaul  policies  on  ship  acquisition 
requirements  as  well  as  on  total  expected  fleet  cost.   Fur- 
ther, we  have  seen  the  importance  of  considering  the  real  cost 
differential  between  scheduled  and  unscheduled  overhauls. 

The  significant  conclusions  to  be  drawn  from  this  thesis 
are  threefold: 

(1)   That  the  use  of  transient  Markov  processes  to  model 
ship  deterioration  and  failure  behavior  is  reasonable. 
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(2)  That  such  models  are  potentially  of  significant  use 
to  Navy  (and  other)  decision  makers. 

(3)  That  a  great  deal  of  initial  effort  is  required  in 
the  area  of  data  analysis  to  develop  material  condition 
indices  and  reasonably  accurate  estimates  of  transition  prob- 
abilities . 
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2  2  0. C90000 

1  3  0.909000 

2  3  0.091000 

1  4  0.909100 

2  A  0.090900 

1  5  0.909090 

2  5  0.090910 

1  6  0.909091 

2  6  0. C90909 

1  7  0.818181 

2  7  0.181818 

1  8  0.636363 

2  8  0.163636 

1  9  0.832726 

2  9  0. 167273 

1  10  0.833454 

2  10  0.166545 

1  11  0.833308 

2  11  0.166691 

1  12  0.833337 

2  12  0.  166662 

1  13  0.749997 

2  13  0.250001 

1  14  0.774999 

2  14  0.224999 

1  15  0.767499 

2  15  0.232500 

1  16  0.769749 

2  16  0.230250 

1  17  0.769074 

2  17  0.230925 

1  18  0.769276 

2  18  0.230722 

1  19  0.692288 

2  19  0.307710 

1  20  0.723083 

2  20  0.276915 

1  21  0.710765 

2  21  0.289233 

1  22  0.715692 

2  22  0.284306 

1  23  0.713721 

2  23  0.286277 

1  24  0.714509 

2  24  0.285488 

1  25  0.642743 

2  25  0.357255 

1  26  0.678626 

2  26  0.321371 

1  27  0.660684 

2  27  0.339313 


1  28  0.669655 

2  28  0.  330342 

1  29  0,665170 

2  29  0.334828 

1  30  0.667412 

2  30  0.332585 
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0.022003 
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25 
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10 

0.000698 

1 

25 

11 

0.000011 
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25 

11 

0.000048 
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25 

12 

0.000000 
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12 
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1 

25 
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0.000000 
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13 

0.000000 

1 
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14 

o.ocoooo 
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14 

0.000000 

1 

25 
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15 
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1 

25 

16 

0.000000 

2 

25 

16 

0.000000 

1 

25 

17 

0.000000 
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25 

17 

0.000000 

I 

25 

18 
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25 
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25 
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25 

20 

0.000000 

2 

25 

20 

o.ocoooo 

1 

25 
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25 
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23 
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1 

25 

25 

0.000000 
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25 

25 
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1 

26 

0 
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2 

26 
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26 

1 

0.002933 
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26 
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0.000191 
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26 

2 

0.0  18635 
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26 
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0.002550 
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26 

3 

0.0*4557 
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26 

3 
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5 
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0.154105 
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26 

6 

0.085806 
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1 

26 

7 

0.085641 
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26 

7 

0.062313 

1 
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8 

0.030342 

2 

26 

8 

0.028872 

1 
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9 

0.006615 
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26 

9 

0.0C8339 
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26 

10 
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0.001429 
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2 
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2 

27 
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2 
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9 
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1 
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10 
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0.000030 
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13 
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16 
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28 

19 

0.000000 

1 

28 

20 

0. ocoooo 

2 
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11 
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11 
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0.000247 
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29 

13 
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29 

13 
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1 

29 

14 
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29 

14 
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29 

15 
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15 
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0.000000 

1 

29 

29 

o.ocoooo 

2 
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0.000000 
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30 
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30 
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30 
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30 
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30 
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1 

30 
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30 

7 
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30 
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30 
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1 

30 
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2 

30 
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1 

30 

10 
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30 
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2 

30 

11 
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30 
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2 

30 

12 
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30 
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30 
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o.ocoooo 
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30 
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30 

19 
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1 

30 
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30 

20 

0.000000 
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30 

21 
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30 

21 
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30 

22 
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30 

22 
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30 

23 
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2 

30 

23 

0.000000 

1 

30 

24 
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2 

30 

24 
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1 

30 

25 

0.000000 

2 

30 

25 

0.000000 

1 

30 

26 

0.000000 

2 

30 

26 

o.ocoooo 

1 

30 

27 

o.ocoooo 

2 

30 

27 

0.000000 

1 

30 

28 

o.ocoooo 

2 

30 

28 

0.000000 

1 

30 

29 

0.000000 

2 

30 

29 

0.000000 

1 

30 

30 

0.000000 

2 

30 

30 

o.ocoooo 
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The  expected  number  of  overhauls 
the  ship  will  experience  is: 

6.37877 
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APPENDIX  B 

Computer  Program  for  Determining 
Required  Fleet  Size, 
Homogeneous  Fleet 
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APPENDIX  C 

Recursion  Relationships  for  Age  and 
Time-Since-Overhaul 
Dependent  Models 
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APPENDIX  C 
RECURSION  RELATIONSHIPS  FOR  AGE  AND 
TIME-SINCE -OVERHAUL  DEPENDENT 
MODELS 

The  following  are  the  recursion  formulas  developed  for 
the  age  and  time-since-overhaul  models  for  both  2  and  4  con- 
ditions : 

1.   The  2  Condition  Model: 

P(k,t,j|k0,to,Jo)  =  P(l,t-l,j-l|ko,to,Jo)p(k|l,t-l,j-l) 

3    =  "f • • » ft 
k  =  1,2 

Since  being  in  condition  2  is  tantamount  to  undergoing  an 
overhaul  which  will  return  the  ship  to  condition  1  with  cer- 
tainty, and  since  j  takes  a  value  of  1  only  upon  the  ship's 
coming  out  of  overhaul, 

ft-1 
P(k,t,j|ko,to,Jo)  =   .JiP(2,t-l,£|k0,to,Jo)   j=l,  k-1 

k  0  j=l,  k=2. 

Also , 

TT(k,t,j,X>  =  7T(l,t-l,j-l,x)p(k|l,t-l,j-l) 
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iT(k,t,j,x)  =  *(i,t-l,j-l)p(k|lft-l,j-l), 


Again, 


7r(k,t,j,x)  =  ■ 


t-1 

I    Tr(2,t-l,A,x)    j=l,  k=l 
1=1 

0  j=l,  k=2 


2.   The  multi-condition  (4)  model,  following  the  same  approach: 

3 

I   P(i,t-l,j-l|k  ,t  ,j  )p(k|i,t-l,j-l) 
i=l  °  o 

t-1 
P(k,t,j|ko,to,Jo)=.|. J  P(4/t-l,£|ko,to,Jo)   j-lf  k=l 

0  j=l;  k=2,3,4. 


Tr(k,t,j,x)  =  d 


I    7r(i,t-l„  j-l,x)p(k|i,  t-1,  j-1)   j=2,...,t; 
i=l  k=l,2,3 


I   7T  (i,t-l,  j-l,x-l)p(k|i,t-l,  j-1)  j=2,...,t; 
i=l  k=4 

t-1 

I    Tr(4,t-l,£,x)    j=l,  k=l 
1=1 

0  j=l;  k=2,3,4. 
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APPENDIX  D 

Determination  of  Transition  Probabilities 
for  the  Example  for  the  Multi- 
Condition,  Age  and  Time- 
Since-Overhaul  Dependent 
Model 
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APPENDIX  D 
DETERMINATION  OF  TRANSITION  PROBABILITIES 
FOR  THE  MULTI-CONDITION,  AGE  AND  TIME- 
SINCE  -OVERHAUL  DEPENDENT  MODEL 

The  available  data  provided  the  dates  of  overhauls  and 
restricted  availabilities  for  CRUDESLANT  ships  from  1963  to 
1971  (approximately) .   Thirty-three  ships  were  selected  for 
which  such  data  was  available  for  a  complete  cycle  (i.e., 
overhaul  to  overhaul) .   These  ships  were  grouped  into  the 
three  age  categories  described  in  Chapter  III.   For  each  cate- 
gory a  histogram  (Figure  D-l)  was  developed  showing  the  total 
number  of  restricted  availabilities  occurring  during  each  time 
unit  after  overhaul  (up  to  three  years  for  category  I  and  four 
years  for  categories  II  and  III  -  since  COMCRUDESLANT  has 
never  operated  with  an  overhaul  policy  with  cycles  in  excess 
of  four  years,  little  data  was  available  for  ships  out  of 
overhaul  for  longer  than  that  period) . 

The  results  of  this  were  then  examined  for  some  logical 
grouping  of  ships  based  on  time  since  last  overhaul.   It 
appears  that  in  each  case  the  number  of  RAV's  during  the  first 
unit  of  time  is  markedly  higher  than  that  for  subsequent  time 
periods.   Further,  there  also  appears  to  be  a  relatively  well 
defined  difference  between  the  second  and  third  periods  and 
later  periods.   Therefore,  the  first  category  of  time  out  of 
overhaul  consists  simply  of  j  =  1,  the  second  category  only  of 
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Category  I 
Ship  Age  1-3  Years,  Inclusive 
11  Ships 
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Category  II 
Ship  Age  4-13  Years,  Inclusive 
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Category   III 
Ship  Age   14-25   Years, Incl. 
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23456789 
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FIGURE  D-l 
Histogram  of  RAV  Frequency 
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j  =  2,3,  while  the  third  therefore  consists  of  j  =  4,...,t 
(unless  j  is  otherwise  restricted  by  some  overhaul  policy) . 
With  time  since  last  overhaul  categorized  thusly,  some 
estimate  of  state  probabilities  was  required.   This  was  done 
by  first  making  the  arbitrary  assumption  that  one  out  of 
every  three  RAV's  (which,  in  general,  represent  unexpected 
equipment  failures  of  a  relatively  critical  nature)  was  due 
to  a  ship's  being  in  condition  4.   Therefore,  one-third  the 
number  of  RAV's  averaged  over  time  and  number  of  ships  (for 
each  category  -  there  now  being  3  x  3,  or  9,  categories), 
provided  an  estimator  of  the  probability  of  being  in  condition 
4.   It  was  then  arbitrarily  assumed  that  the  state  probabili- 
ties for  condition  3  were  equal  to  those  for  condition  4,  and 
that  those  for  condition  1  equaled  those  for  condition  2. 
Since  these  must  sum  to  unity,  these  assumptions  define  the 
state  probabilities.   Working  "backwards"  from  the  estimates 
of  state  probabilities,  estimators  for  transition  probabili- 
ties were  developed,  with  the  primary  assumption  being  that 
a  ship  would  be  more  likely  to  remain  in  the  same  condition 
over  a  transition  than  it  would  be  to  change  to  any  other 
given  condition.   The  matrices  of  these  transition  probability 
estimators  comprise  Table  (4)  in  Chapter  III. 

The  high  "failure  rate"  immediately  following  overhaul 
is  somewhat  counterintuitive.   It  is  felt  by  the  author  that 
this  is  due  primarily  to  two  factors.   First,  when  a  ship  is 
overhauled,  disruption  to  major  mechanical  or  electrical/ 
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electronic  systems  is  frequently  involved.   Such  disruptions 
may  in  some  cases  result  in  casualties  to  these  systems  once 
operation  is  resumed.   Secondly,  what  is  administratively 
described  as  a  restricted  availability,  immediately  following 
overhaul,  may  in  fact  be  a  continuation  of  the  overhaul  to 
complete  one  or  more  unfinished  jobs.   This  could  result  from 
official  disfavor  toward  failing  to  complete  an  overhaul  on 
schedule. 

Certain  of  the  assumptions  made  in  developing  the  transi- 
tion probability  matrices  would  not  be  tenable  in  actually 
applying  the  model.   However,  for  the  results  of  the  model  to 
be  reasonably  valid,  it  is  important  that  the  initial  values 
estimated  for  transition  probabilities  be  based  on  complete, 
accurate  data  insofar  as  possible,  and  the  criteria  for  esti- 
mating a  ship's  material  condition  be  applied  in  a  consistent 
manner  to  all  ships.   As  noted  in  Chapter  II,  while  efforts 
continue  to  develop  a  material  condition  index  which  will 
permit  a  quantitative  description  of  a  ship's  material  condi- 
tion, no  such  usable  index  currently  exists.   Some  reasonable 
approximation  of  such  an  index  is  required,  then.   Of  signifi- 
cant value  in  developing  such  an  approximation  might  be  the 
results  of  type  commanders'  material  inspections,  or  inspec- 
tions by  the  Board  of  Inspection  and  Survey  (INSURV  Board) . 
Other  possible  parameters  that  could  be  used  in  estimating  a 
ship's  material  condition  include  RAV  frequency  (as  in  the 
example),  CASREPT  data  (frequency,  severity,  etc.),  and  the 
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number  of  man-hours  expended  in  maintenance  or  the  volume  of 
deferred  maintenance  (under  the  "3-M"  system)  reported  by  a 
ship.   Here  one  should  be  careful,  however,  as  inconsistencies 
occur  from  ship  to  ship  in  the  accuracy  and  thoroughness  of 
the  reporting  of  maintenance  data. 

Whatever  combination  of  parameters  is  decided  upon,  the 
best  approach  to  estimating  transition  probabilites  is  prob- 
ably to  choose  several  representative  ships  and  use  past  data 
to  "track"  them  over  time,,  observing  how  they  changed  in 
material  condition  and  noting  age  and  time  since  overhaul  at 
the  time  changes  took  place. 

It  should  be  obvious  that  once  transition  probability 
estimators  have  been  determined  and  a  model  such  as  this 
implemented,  the  estimators  (as  well  as  the  criteria  for 
grouping  ships  into  categories)  should  be  updated  as  behavior 
of  ship  material  condition  over  time  is  actually  observed. 
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APPENDIX  E 

Computer  Program  for  the  Multi- 
Condition,  Age  and  Time-Since 
Overhaul  Dependent  Model 
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APPENDIX  F 

Aggregated  State  Probabilities  ($(k,t)) 
for  the  Multi-Condition,  Age  and 
Time-Since-Overhaul  Dependent 
Model  Example 
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k  t  4»(k,t) 

1  L  1.000000 

2  1  0,000000 

3  1  C. 000000 

4  1  0.000000 

1  2  0.505000 

2  2  C. 253000 

3  2  C. 121000 

4  2  C. 121000 

1  3  C. 554177 

2  3  C. 332377 

3  3  0.080923 

4  3  C. 032523 

1  4  0.489522 

2  4  0.386187 

3  4  0.078330 

4  4  0.045961 

1  5  0.503031 

2  5  C. 420550 

3  5  0.053875 

4  5  0.022543 

1  6  0.488774 

2  6  C. 440040 

3  6  C. 046367 
A  6  C. 024817 

1  7  C. 490047 

2  7  C. 449073 

3  7  0.039713 

4  7  C. 021166 

1  6  0.487707 

2  8  0.453824 

3  8  0.037176 

4  8  0.021291 

1  9  P. 487845 

2  9  C. 456134 

3  9  0.035444 

4  9  C. 020574 

1  10  0.402857 

2  10  C. 374581 

3  U  C. 118369 

4  10  C. 104191 

1  11  0.436053 

2  11  C. 325366 

3  11  0.142963 

4  11  C. 095615 

1  12  0.420386 

2  12  0.304877 

3  12  C. 165959 

4  12  C. 103774 
1  13  0.411193 
?  13  C. 298589 

3  13  0.168298 

4  13  C. 101915 

1  14  0.430217 

2  14  0.300714 
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3  14  0.168191 

4  14  0.100872 

1  15  C. 431027 

2  15  0.302202 

3  15  0.167020 

4  15  0.099744 

1  16  0.431010 

2  16  0.303340 

3  16  0.166226 

4  16  0.099418 

1  17  C. 430644 

2  17  0.303499 

3  17  0.166170 

4  17  0.099680 

1  18  0.430686 

2  18  0.303412 

3  18  0.166181 

4  18  0.099713 

1  19  C. 430589 

2  19  C. 303291 
?  19  C. 166292 
4  19  C. 099819 

1  20  0.430624 

2  2C  0.303219 

3  20  C. 166332 

4  20  0.C99815 

1  21  C. 430622 

2  21  C. 303202 

3  21  0.166349 

4  21  0.099816 

1  22  C. 430629 

2  22  0.303202 

3  22  0.166348 

4  22  C.C99809 

1  23  0.430631 

2  23  0.303208 

3  23  0.166344 

4  23  C. 099804 

1  24  C. 430630 

2  24  C. 303212 

3  24  0.166341 

4  24  C. 099804 

1  25  C. 430630 

2  25  0.303213 

3  25  0.166340 

4  25  0.099803 

1  26  C. 430628 

2  26  C. 303212 

3  26  0.166340 

4  26  0.099804 

1  27  C. 430628 

2  21  0.303212 

3  27  0.166340 

4  27  C. 099804 
1  28  C. 410627 
?  28  0.303211 

3  28  0.166340 

4  28  0.099804 

1  29  C. 430627 

2  29  0.303211 
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3  29  C. 1663*0 

4  29  0,09980* 

1  30  0.430627 

2  30  0.303211 

3  30  0.1663*0 

4  30  0.099804 

1  31  C. 430626 

2  31  0.303210 

3  31  0.166340 

4  31  0.099803 

1  32  C. 430625 

2  32  0.303210 

3  32  C. 166339 

4  32  0.C99803 

1  33  C. 430625 

2  33  0.303210 

3  33  0.166339 

4  33  0.0<>9803 

1  34  0.430624 

2  34  C. 303209 

3  34  C. 166339 

4  34  0.099803 

1  35  0.430624 

2  35  0.303209 

3  35  0.166339 

4  35  0.C99803 

1  36  0.430623 

2  36  0.303203 

3  36  0.16633R 

4  36  0.09<?803 

1  37  0.430623 

2  37  0.303208 

3  37  C. 166338 

4  37  0.099802 

1  38  C. 430622 

2  38  0.303208 

3  38  C. 166338 

4  38  C. 099802 

1  39  0.430621 

2  39  0.303207 

3  39  0.166338 

4  39  0.099802 

1  4G  0.476233 

2  40  0.336843 

3  4C  0.126713 

4  4C  0.060179 

1  *l  0.464512 

2  41  0.359954 

3  41  0.113093 

4  41  0.062408 

1  42  0.464399 

2  42  0.367C28 

3  42  0.106888 

4  42  0.061651 

1  43  0.462283 

2  43  0.368798 

3  43  0.105820 

4  43  0.063065 

1  44  0.461337 

2  44  0.367906 
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3  44  0.106525 

4  44  0.064197 

1  45  0.461803 

2  45  0.367422 

3  45  0.106675 

4  45  0.064065 

1  46  0.461846 

2  46  0.367347 

3  46  0.106720 

4  46  0.064050 

1  47  0.461945 

2  47  0.367399 

3  47  0.106654 

4  47  0.063966 

1  48  0.461^53 

2  48  0.367470 

3  46  0.106600 

4  48  0.063939 

1  49  0.461935 

2  49  0.367495 

3  49  0.106585 

4  49  0.063946 

1  5C  0.461931 

2  5C  0.367499 

3  50  0.106583 

4  5&  0.063949 

1  51  0.461927 

2  51  0.367496 

3  51  0.106585 

4  51  C. 063952 

1  52  0.461927 

2  52  0.367493 

3  52  0.106586 

4  52  0.063952 

1  53  C. 461928 

2  53  0.367492 

3  53  0.106586 

4  53  0.063952 

1  54  0.461927 

2  54  0.367492 

3  54  0.106586 

4  54  0.063952 

1  55  C. 461927 

2  55  0.367492 

3  55  0.106586 

4  55  0.063952 

1  56  0.461927 

2  56  0.367492 

3  56  C. 106586 

4  56  0.063951 

1  57  0.461926 

2  57  0.367491 

3  57  0.106585 

4  57  C. 063951 

1  58  0.461925 

2  58  0.367490 

3  58  0.106585 

4  58  0.063951 

1  59  0.«»61925 

2  59  0.367491 
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3  59  0.106585 

4  59  C. 063951 

1  60  0.461925 

2  60  0.367490 

3  60  0.106585 

4  60  C.063951 

1  61  0.461924 

2  61  0.367490 

3  61  0.106585 

4  61  0.063951 

1  62  C. 461924 

2  .  62  C. 367489 

3  62  0.106584 

4  62  C.063951 

1  63  0.461923 

2  63  0.367489 

3  63  0.106585 

4  63  C.063951 

1  64  0.461922 

2  64  0.367488 

3  64  0.106584 

4  64  0.063951 

1  65  0.461922 

2  65  C. 367488 

3  65  0.106584 

4  65  0.063951 

1  66  C. 461921 

2  66  0.367487 

3  56  0.106584 

4  66  0. 063951 

1  67  0.461921 

2  67  0.367487 

3  67  C. 106584 

4  67  0.063951 

1  68  0.461920 

2  68  C.367486 

3  68  0.106584 

4  68  0.063951 

1  69  0.461919 

2  69  0.367486 

3  69  C. 106584 

4  69  C. 063950 

1  7C  0.461919 

2  7C  C.367486 

3  7C  0.106583 

4  7C  0.063950 

1  71  0.461918 

2  71  0.367485 

3  71  C. 106583 

4  71  0.063950 

1  72  C. 461918 

2  72  C. 367485 

3  72  C. 106583 

4  72  C.063<»50 

1  73  0.461917 

2  73  0.367484 

3  73  0.106583 

4  73  0.063950 

1  74  0.461917 

2  74  0.367484 
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74 
74 
75 
75 
75 
75 


0.106583 
0.063950 
0.461916 
0.367483 
0.106582 
0.063950 
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APPENDIX  G 

Computer  Program  for  Determining 
Required  Fleet  Size,  Mixed  Fleet 


137 


******* 


< 

CO 

z 


u. 

o 

co 

z 
o 


<  < 

—    CC 

co  o 

a.  z 

a 

o  UJ 

x 


o 

UJ 

CO 

HH 

(X 

u 
co 
ai 
O 

co 


co 

O 

•— ■ 
ct 

< 

;> 


h-  at. 


o 


oo 


cl    ►►- 

O  h-  •— 
UJ  2T 
</)  UJ  •-"• 
UJ  _J  _l 
i-t  u_ 


•  < 

Q  cc 

x  z 

•"•  !JJ 

5:  o 


aj 

CD 

o 

CC 

Cl 

co 

UJ 


UJ 
UJ  X 
X  K 


cc  o 
o 

u.  o 

UJ 

o  O-  < 

•—  cr 


3 
a 


o 

a 

CO  Z 

•— • 

3 

CD 

h- 

o 

-J    "H 

a 

UJ 

<   h- 

O 

o 

Z  < 

o 

O  31 

o 

a' 

•— •  >— « 

o 

CL 

K  X 

a: 

«i  a 

a 

CO 

CC   (X 

_4 

uj  a 

«• 

X 

a.  a 

»- 

K 

o  <r 

_j 

u_ 

X.   * 

* 

*  * 

ar  "v 

V. 

>»  *v 

CO 

uj  a. 

x  ♦-« 

H-  X 

to 

UJ 

<l  < 

3 
H-  O 
3  •-• 
Cl  > 

Z  ►-« 

M    Q 

z 

CO  >— 

< 

Q  O 
UJ  U- 

CC 

»~"  -J 
3  Z> 

o< 

UJ  X 

«:  cc 

UJ 

> 
•  o 
> 
«  z 

•— ■ 

a 

UJ  o 
H-  Z 

<  UJ 
X  CO 

o 

UL 

U-  O 

a 
to 

O  UJ 


o  •-• 


UJ  < 
CO  CO 

o 

Z  CC 

»-.  a. 


< 


o 

X 


<l 
a 

-j 

LL 

3 


O 
UJ 
X 

•— I 
U- 


co 

z 

z 

CO 


CM 


—  O 

CO  UJ 

<N  X 

i— I  *-* 

cr. 

CL  2 


CJ  *-* 

U   CO 


to 

OJ 

r-t 

o 


II 

z 

O  cr> 

*~  V 

Z  CO 

CC  o 

Q-  O 


to 

>- 

CO 


CO 

CC 

a 

CC 


CO 
CC 

a 

•> 
oj 

CC 

a 


# 

CC 
CL 

* 

z 
+ 

co 
cc 
a 

I 
o 


* 

co 
cc 

UL 

* 

m 
z 
•»  + 
in  <— 
cc  f\J 
a.  cc 

*    CL 
CO     | 

z  o 

♦  • 

cc  — 

a  * 

*  rvj 
4-  cc 

Z  O. 

♦  * 
m  c\j 

cc  z 

a  + 

#  — 

co  i-i 


o  -* 

CM    + 


f\)  cO  «j-  tr\  cc 


♦  +  + 


O  O   +   4 


CO 

»-  O 
_)    II 


UJ  u 


•f  ♦  ♦ 

_    M     ►_ t    <J 


O  cc 

ii  a. 

"5  II 

a  cc 

O  CL 


>o 
a 

o 

••  H 

CL 
w    w    w    w    »_    »-•    -4" 

cc  cc  cc  cc  <  *    II 

CL    CL   CL    OL  O   CO  r-i 
II      II     II      II  Z 

<M  m  ^  u\  i-  i- 

aaaajDu 

tx    a  CL    U.  CL    CL  Q 


O  vO  >o  >o 

DUOC 


't  -4-  «4"  < 

II  II  II  || 

rvj  ro  •>*•  crv 

z  z  z  z 

uuuu 

DUQU 


CO 

cc 
a. 

I 

o 

-•  o 

^-  OJ 

*   I 

CO  CO 

cc  z 

CL    •¥ 
*    ,* 

in  Z 
Z  + 

+  ro 


< 
> 


CC   -» 

o  o 
•^  a 


co 
z 

« 

z 
m 


r\j 


Z  « 


LJ 

z 

IJJ 


o 
z 

co 


Q 

Z 
UJ 


•^  Z 

Z    II 

ii  or 
3  <r 
ac  > 


rsj 

^. 

■¥ 
—I 

z 
II 

CO 


CC 


<r  <i  cl  o 

K  I-  •-  Z 

<  O  ^  UJ 

o  a  co 


_j  a. 
<i  "i 


3  3  J 
CL   CL    CL 


c 

Z 

cu 


*  *  * 


138 


— 

o 

ITt 

-H 

•  »■ 

~- 

«_ 

h- 

O. 

a 

_J 

►— < 

M-4 

LL 

* 

* 

X. 

00 

•  » 

to 

•  » 

i 

t- 

n 

t- 

O  D 

3  ^ 

Z> 

^ 

^ 

u. 

U' 

a 

LU 

u. 

Thesis 

135357 

E3125 

Eldred 

Development 

of 

a  com- 

parat 

ive   Markov  model    for 

ship 

overhaul 

po 

1  icies. 

thesE3125 

°um,  °Pmem  °' 3  comPara"«  Markov  mode 


3  2768  001  89271  4 

DUDLEY  KNOX  LIBRARY 


